Google 



This is a digital copy of a book thaï was prcscrvod for générations on library shelves before it was carefully scanned by Google as part of a project 

to make the world's bocks discoverablc online. 

It has survived long enough for the copyright to expire and the book to enter the public domain. A public domain book is one that was never subject 

to copyright or whose légal copyright term has expired. Whether a book is in the public domain may vary country to country. Public domain books 

are our gateways to the past, representing a wealth of history, culture and knowledge that's often difficult to discover. 

Marks, notations and other maiginalia présent in the original volume will appear in this file - a reminder of this book's long journcy from the 

publisher to a library and finally to you. 

Usage guidelines 

Google is proud to partner with libraries to digitize public domain materials and make them widely accessible. Public domain books belong to the 
public and we are merely their custodians. Nevertheless, this work is expensive, so in order to keep providing this resource, we hâve taken steps to 
prcvcnt abuse by commercial parties, including placing lechnical restrictions on automated querying. 
We also ask that you: 

+ Make non-commercial use of the files We designed Google Book Search for use by individuals, and we request that you use thèse files for 
Personal, non-commercial purposes. 

+ Refrain fivm automated querying Do nol send automated queries of any sort to Google's System: If you are conducting research on machine 
translation, optical character récognition or other areas where access to a laige amount of text is helpful, please contact us. We encourage the 
use of public domain materials for thèse purposes and may be able to help. 

+ Maintain attributionTht GoogX'S "watermark" you see on each file is essential for informingpcoplcabout this project and helping them find 
additional materials through Google Book Search. Please do not remove it. 

+ Keep it légal Whatever your use, remember that you are lesponsible for ensuring that what you are doing is légal. Do not assume that just 
because we believe a book is in the public domain for users in the United States, that the work is also in the public domain for users in other 
countiies. Whether a book is still in copyright varies from country to country, and we can'l offer guidance on whether any spécifie use of 
any spécifie book is allowed. Please do not assume that a book's appearance in Google Book Search means it can be used in any manner 
anywhere in the world. Copyright infringement liabili^ can be quite severe. 

About Google Book Search 

Google's mission is to organize the world's information and to make it universally accessible and useful. Google Book Search helps rcaders 
discover the world's books while helping authors and publishers reach new audiences. You can search through the full icxi of ihis book on the web 

at |http: //books. google .com/l 



Google 



A propos de ce livre 

Ceci est une copie numérique d'un ouvrage conservé depuis des générations dans les rayonnages d'une bibliothèque avant d'être numérisé avec 

précaution par Google dans le cadre d'un projet visant à permettre aux internautes de découvrir l'ensemble du patrimoine littéraire mondial en 

ligne. 

Ce livre étant relativement ancien, il n'est plus protégé par la loi sur les droits d'auteur et appartient à présent au domaine public. L'expression 

"appartenir au domaine public" signifie que le livre en question n'a jamais été soumis aux droits d'auteur ou que ses droits légaux sont arrivés à 

expiration. Les conditions requises pour qu'un livre tombe dans le domaine public peuvent varier d'un pays à l'autre. Les livres libres de droit sont 

autant de liens avec le passé. Ils sont les témoins de la richesse de notre histoire, de notre patrimoine culturel et de la connaissance humaine et sont 

trop souvent difficilement accessibles au public. 

Les notes de bas de page et autres annotations en maige du texte présentes dans le volume original sont reprises dans ce fichier, comme un souvenir 

du long chemin parcouru par l'ouvrage depuis la maison d'édition en passant par la bibliothèque pour finalement se retrouver entre vos mains. 

Consignes d'utilisation 

Google est fier de travailler en partenariat avec des bibliothèques à la numérisation des ouvrages apparienani au domaine public et de les rendre 
ainsi accessibles à tous. Ces livres sont en effet la propriété de tous et de toutes et nous sommes tout simplement les gardiens de ce patrimoine. 
Il s'agit toutefois d'un projet coûteux. Par conséquent et en vue de poursuivre la diffusion de ces ressources inépuisables, nous avons pris les 
dispositions nécessaires afin de prévenir les éventuels abus auxquels pourraient se livrer des sites marchands tiers, notamment en instaurant des 
contraintes techniques relatives aux requêtes automatisées. 
Nous vous demandons également de: 

+ Ne pas utiliser les fichiers à des fins commerciales Nous avons conçu le programme Google Recherche de Livres à l'usage des particuliers. 
Nous vous demandons donc d'utiliser uniquement ces fichiers à des fins personnelles. Ils ne sauraient en effet être employés dans un 
quelconque but commercial. 

+ Ne pas procéder à des requêtes automatisées N'envoyez aucune requête automatisée quelle qu'elle soit au système Google. Si vous effectuez 
des recherches concernant les logiciels de traduction, la reconnaissance optique de caractères ou tout autre domaine nécessitant de disposer 
d'importantes quantités de texte, n'hésitez pas à nous contacter Nous encourageons pour la réalisation de ce type de travaux l'utilisation des 
ouvrages et documents appartenant au domaine public et serions heureux de vous être utile. 

+ Ne pas supprimer l'attribution Le filigrane Google contenu dans chaque fichier est indispensable pour informer les internautes de notre projet 
et leur permettre d'accéder à davantage de documents par l'intermédiaire du Programme Google Recherche de Livres. Ne le supprimez en 
aucun cas. 

+ Rester dans la légalité Quelle que soit l'utilisation que vous comptez faire des fichiers, n'oubliez pas qu'il est de votre responsabilité de 
veiller à respecter la loi. Si un ouvrage appartient au domaine public américain, n'en déduisez pas pour autant qu'il en va de même dans 
les autres pays. La durée légale des droits d'auteur d'un livre varie d'un pays à l'autre. Nous ne sommes donc pas en mesure de répertorier 
les ouvrages dont l'utilisation est autorisée et ceux dont elle ne l'est pas. Ne croyez pas que le simple fait d'afficher un livre sur Google 
Recherche de Livres signifie que celui-ci peut être utilisé de quelque façon que ce soit dans le monde entier. La condamnation à laquelle vous 
vous exposeriez en cas de violation des droits d'auteur peut être sévère. 

A propos du service Google Recherche de Livres 

En favorisant la recherche et l'accès à un nombre croissant de livres disponibles dans de nombreuses langues, dont le français, Google souhaite 
contribuer à promouvoir la diversité culturelle grâce à Google Recherche de Livres. En effet, le Programme Google Recherche de Livres permet 
aux internautes de découvrir le patrimoine littéraire mondial, tout en aidant les auteurs et les éditeurs à élargir leur public. Vous pouvez effectuer 
des recherches en ligne dans le texte intégral de cet ouvrage à l'adresse fhttp: //book s .google . coïrïl 




,<'^:; 



THEEK GENT 






UNIVERSITEITSBIB 




9000000 




I 
I 



■ IIPIH. Ç 



\ 



NOUVELLES METHODES 

POUR LA DÉTERMINATION / 

DES f4^|//0^ 

ORBITES DES COMÈTES; 



PAR A. M. LEGENDRE, 

Membre de Plnstitut et de la Légion d'honneur , de la Société 

rojale de Londres , &c. 




A PARIS, 

Chez FiRMiN DIDOT, Libraire pour les Mathématiques , la Marine , 
rArchitecture , et les Éditions stéréotypes , rue de Thionville, n" n 6. 

AN XI II — i8o5. 



'/■ -■■■ ' 

I 



..''' 



f 4 "l. 



NOUVELLES MÉTHODES 

I 

POUR tA DÉTERMINATION 

DES ORBITES DES COMÈTES. 



U 



£jE problème dont je m^occuperai dans ce Mémoire , consiste à 
'déterminer Porbite d'ane comète d'après trois observations don- 
nées de sa longitude et de sa latitnde. Depuis Newton , qui le 
premier a donné pour sa solution des constructions géométriques 
fort ingénieuses , ce problème a fait successivement Fobjet des 
recherches d'an grand nombre de géomètres. 

Parmi ceux qui Tont traité avec le plus de succès^ on doit 
distinguer particulièrement Lambert, qui a donné de très-beaux 
théorèmes sur le mouvement des planètes et des comètes daus 
son ouvrage intitulé : Insigniores orbitœ cometarum proprieêates, 
et dans les Mémoires de PAcàdémie de Berlin , année 1771. 

Dans ce même recueil y année 1778 , La Grange a discuté les 
principales méthodes connues jusqu'alors , et après avoir- fait 
connoître les causes de leur imperfection , il a donné l'analj'se 
complète du problème^ fondée sur une belle théorie a laquelle 
il a ajouté ensuite de nouveau^: développemens dans le volume 
de 1783. Cette méthode n'auroit sans doute rien laissé à désirer, 
si son illustre auteur exi eAt fait V^ppHoation à des exemples , 
ce qui l'auroit conduit lui*mème h y apporter les modifications 
nécessaires pour en rendre Tusaga ftioile dans la pratique. 

La Place , dans les Hémoiri?é de TAcadémie des 3cienoea de 
Paris , année 1780 , a proposé une autre méthode qui revient ji 
supposer connue une portion infiniment petite de la trajeetoire 
apparente de la comète « <$t à conclure de cette portion la gran- 
deur et la position de la trajectoire vraie. Far la supposition 
des quantités infiniment pQtite^ y les formules se simplifient et 
conduisent à une solution peu compliquée; mais la difficulté est 
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de déterminer avec précision les coëfficiens différentiels du pre- 
mier et du second ordre , tant de la longitude que de la latitude. 

Il semble , au premier coup-d'œil , qu'en liant ensemble plu* 
sieurs observations par la méthode des interpolations , on en 
déduira les coëfficiens dont il s'agit avec d'autant plus d'exac- 
titude qu'il y a plus d'observations combinées ; et cela auroit 
lieu en effet, si les observations étoient exemptes d'erreur, ou 
si les lieux de la comète étoient calculés d'après une formule 
exacte ; mais il en est autrement dans l'état réel des choses. 
Comme toutes les observations sont affectées de quelqu'erreur , 
et que cette erreur ne suit aucune loi d'une observation à l'autre, 
il s'ensuit que plus on combinera d'observations, et plus l'erreur 
des coëfficiens différentiels qui en sont déduits pourra devenir 
sensible. 

En effet , considérons trois observations de longitude a' y a% a% 
faites à des intervalles de temps que pour plus de simplicité 
nous supposerons égaux à l'unité. Ces longitudes répondront 
aux temps r— 1,0, +1, et la longitude pour un temps quel- 
conque t compris entre — i et + i , aura pour expression 
générale , 

^ = a' + i(a'^— a')/ + i(a'" — 2a'+a')^; 
d'où l'on déduit les coëfficiens différentiels 

dt '^ '' 

Considérons ensuite les cinq longitudes o,_<^f a", a% a"^, qui 
répondent pareillement aux temps -—3,' — i, o, +i, +a, res-. 
pectivement ; on en déduira la longitude au bout du temps /, 
X = a"+At + B <• + C *" + &c. , où l'on aura 

A ou ^ = f («"'-«') - h (cf^-d) 
.3B ou ^ = f (a'"— aa'+a') — T^(a"— 20"+ a). 
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Supposons maintenant que Terreur sur la longitude a soit ^a'i 
alors , suivant les équations (a) , les erreurs des coëffîciens diffé- 
rentiels j daes à cette cause , seront 

Mais par les équations ( 6 ) les erreurs de ces coefficiens seroient 

d'où Ton voit que dans le cas de cinq observations combinées ^ 
Terreur des coëffîciens différentiels , due à la cause mentionnée , . 
est plus grande dans le rapport de 4 à 3 , que celle qui a lieu 
dans le cas de trois observations. Elle augmenteroit encore si 
on combinoit ensemble plus de cinq observations. 

D'après ces réflexions, f ai pensé que ce qu'il y avoit de mieux 
à faire dans le problême des comètes , étoit de partir des données 
immédiates de l'observation , et d'employer tous les moyens 
pour simplifier autant qu'il est possible, les formules et les 
équations qui servent à déterminer les élémens de l'orbite. C'est 
l'objet que je me suis proposé dans ce Mémoire. , 

Je l'ai divisé en deux parties : la première comprend Tanalyse 
générale du problème , avec deux applications détaillées de la 
solution qui en résulte aux comètes de 1781 et de 1769. 

Dans Tanalyse je suppose y oomme cela est indispensable, que 
les trois observations données ne comprennent pas un inter- 
valle de temps de plus de i5 à ao jours , afin que les séries qui 
expriment les coordonnées, tant de Torbite delà comète que de 
celle- de la terre , soient suffisamment convergentes , et qu'on ne 
soit pas obligé d'employer les termes qui contiennent des puis- 
sances du temps supérieures à la troisième. 

Après avoir mis les équations générales du problême sons la 
forme la plus simple dont elles paroissent susceptibles , je déve- 
loppe en particulier, avec beaucoup d'étendue, le cas où les 
trois observations sont faites à des intervalles de temps égaux. 
Cette supposition , qui ne restreint guère la généralité du pro- 
blême, simplifie beaucoup les formules, et contribue même à 
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les rendre plaa exactes , par la disparition cie plusieurs termes 
dont il faudroit tenir compte y si les interralies de temps «itre 
les observations n'étoient pas égaux. 

Lorsqu'on ne fait aucune supposition sur la nature de Torbite, 
le problême offre préqisément autant d'équations que d'incon- 
nues ; mais alors il y a une circonstance dans laquelle les for- 
mules pourroient ne pas donner des résultats assez exacts. C'est 
lorsque l'orbite apparente est située à très- peu près dans le plan 
d'un même grand cercle , et dans ce caa , la distance de la 
comète au soleil , lors de l'observation moyenne y diffère tou-* 
)our& très* peu de la distance de la terre au soleil Si les méthodes 
analytiques ont peu de sucoèa dans ce cas particulier , il en 
résulte au moins une\conttoissance utile sur la distance de la 
comète au soleil , laquelle servira toujours a diriger les premiers 
essais des calculateurs» D'ailleurs , en choisissant trois autres 
observations, à quelque distance des premières , on évitera 
l'inconvénient de ce. cas particulier , à moins qu'il ne se iren^* 
contre tout à* la- fois que la comète smt très-voisine du péri* 
hélie , et que la distance périhélie diffère très*peu de la distance 
de la terre au soleil ; circonstances qui , tout en faisant excep- 
tion, avanceroient beaucoup vers la cannoissance de la véri« 
table orbite. 

Mais comme on suppose communément , d'après les résultats 
de l'observation y que l'orbite est parabolique , cette condition 
donne une équation de plus que d'inconnues, et on a la faculté 
de choiair entre les diverses combinaisons des équations , celle 
qui doit conduire aux résultats les plus exacts. Je suis entré à 
ce sujet dans une discussion fort étendue. J'ai fait voir quelles 
sont les équations qui mèïieroient, dans certains cas, à des 
résultats défectueux , et quelles sont celles sur lesquelles on peut 
établir la solution la moins sujette à être affectée des erreurs 
des observations. 

JOs équations dont il s'agit donnent immédiatement par leur 
résolution les distances de la comète au soleil et à la. terre ; il 
faut ensuite en conclure les élémens de l'orbite. Je donne pour 
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cet éQét toutes les formules uécessaires , et j'ajoute les moyens 
de rfBConnoitre avec certitude si le mouvement est direct out 
rétrograde j si la comète marche vers le périhélie > ou si elle a 
déjà passé par ce point ; si le nœud dont on a calculé la longi- 
tude est le nœud descefndant ou le nosud ascendant ;. de sorte 
qu'avec toutes ces directions les calculs pourront s'exécuter en 
quelque sorte mécaniquement, sans qu'on ait lieu de craindre 
de se tromper sur le signe d'aucune quantité , ou sur la position 
d'aucun des points qu'il importe de déterminer. 

Après avoir traité complètement le cas où les intervalles 
de temps entre les observations sont égaux, il falloit au moins 
tenter de Vaincre les difficultés d'analyse que présente le pro«* 
bléine considéré dans toute Aa généralité. Je suis donc revenu 
sur les équations générales dû ptoMéme , et après leur avoir 
iait subir différentes réductions successives , j'ai trouvé asses 
heureusement qu'elles étcnent sufeqeptibles d'une solution près* 
que aussi simple que celte du premier cas ; de sorte que sans 
aucune interpolation , on pourra appliquer immédiatement le 
calcul à trois observations données , quels que soient les inter^ 
valles de temps qui les séparent 

Tenant ensuite aux applications qui sont relatives à la secondé 
romète de 1781 et à la comète de 176g , j'entre dans les détails 
qui peuvent jeter du jour sur toutes les applications en général* 
J'examine successivement les divers systèmes d'équations qu^n 
peut former y j'en donne la résolution ; je développe les calculs 
principaux avec beaucoup d'étendue, et en général avec plus 
de précision qu'il n'est nécessaire dans les applications ordi^ 
naires:mais celles-ci étant données pour exemple / et devant 
faire juger du degré d'exactitude de la méthode, on ni^pouvoit 
y mettre trop de soin. Les résultats , au surplus , en sont très* 
satisfaisans , et on verra que par cette méthode il est possible de 
déduire de trois observations données d'une comète, dans un 
adsez petit espace de temps, des valeurs fort approchées des 
vrais élémens de son orbite. 

Cependant la connoissance de la véritable orbite ne sera 
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acquise d'une manière certaine , que lorsqu'on aura pu satisfaire 
à trois ou à un plus grand nombre d'observations , dont les 
époques soient assez éloignées entre elles. 

La seconde partie contient les méthodes nécessaires pour 
corriger , à l'aide de ces observations , les élémens connus par 
une première approximation. 

Je prends pour exemple trois observations de la comètp 
de 1769 , faites à de grands intervalles de temps , et je fais voir 
comment y par l'emploi des corrections indéterminées (que 
j'avois déjà indiqué dans les JVfém. de l'Acad: des Se. anu. 1787)9 
on peut parvenir à déterminer la parabole qui satisfait le plus 
exactement possible aux observations. I^a méthode des correc<- 
tions indéterminées a l'avantage de donner , par un seul a^lcul 
et sans aucun tâtonnement , le résultat de toutes les hypothèses 
peu différentes entre elles qu'on pourroit former sur la valeur 
des élémens. Il faut ensuite, lorsque toutes les conditions du 
problême sont exprimées convenablement j déterminer le3 
coëfficiens de manière à rendre les erreurs les plus petites qu'il 
est possible. 

Pour cet effet 9 la méthode qui me paroît la plus simple et 
la plus générale , consiste à rendre minimum la somme des 
quarrés des erreurs. On obtient ainsi autant d'équations qu'il y 
a de coëfficiens inconnus j ce qui achève de déterminer tous les ' 

élémens de l'orbite. 

Comme la méthode dont je viens de parler , et que j'appelle ; 

Méthode des moindres quarrés, peut être d'une grande utilité j 

dans toutes les questions de physique et d'astronomie où il j 

s'agit de tirer de l'observation les résultats les plus exacts qu'elle 
peut offrir ; j'ai ajouté , dans une appendice, des détails parti- 
culiers sur cette méthode , et j'en ai donné l'application à la 
mesure de la méridienne *de France , ce qui pourra servir de 
complément à ce que j'ai déjà publié sur cette matière. 
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analyse du problême. 
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L Soit SAT ( fig. i) le plan dd récliptique , S le soleil , AT 
l'orbité'de lar terre y SA le rayon dirigé vers le premier pokit- 
d^aries supposé' fixe; soient T et C les lieux de la terœ et de- 
là comète , en un même temps t. Du point C abaissez CK per- 
pendiculaire sur le plan de l'écliptique j menez KF, TË/ 
perpendiculaires , et TG parallèle à SA; enfin joignez SC , TC. 

Nous appellerons a?, ^, i , les trois Coordonnées SF, FK, KC 
de la comète, p son rayon Tecteur SC , f sa distance à ht terra 

TC; '^ - 

X , Y, les deux coordonnées S E , E T de la terre , Y son 
rayon vecteur S T ; 

« et ^ la longitude et la latitude géocentrique de la comète y 
c'est-à-dire , les angles GTK , CTK. 
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n. Si Ton prend pour unité la somme des masses du soleil 
et de la terre, et qu'on néglige par rapport à cette somme, la 
différence , d'ailleurs inconnue , £nt|*e la, mas^e de la comète 
et celle de la terre; les* ^équationà différentielles du mouvement 
de la comète seront: 



ddx X ddy ' y ddz' 



z 



dt^ p' ' de: . y' ' d^ r" i/"' 

et celles du mouvement de la terre y (i) 

' 'iàx^ • X ddY ^ Y • 



de V ' dt 



3 > ^*a "^ fr3' 



Ces équations supposent en outre que le tçmps t est mesuré par 
le moyen mouvement du soleil réduit en parties du rayon y 




(a) 

et qu'on fait abstraction de l'action réciproque entre la terre 
et la comète. 

m. Cela posé , soient encore m^ n , p , r, les valears respec« 
tives de s, y, z, v/ et M , N, R, eettes^ de X, Y, V, lorsque 
t = o. Si on suppose qu'à compter de l'époque où < = o , les 
intervalles, dé tempa ne sotent pas troip aonsidémbles en deçà 
el aurdelà de celte époque ;; lea caordunj\âea, taAt de la ccinète 
que ée la barre , seront déterminées arti». uacr axacUitade suffir» 
•aHta par ka équationa aoÎTantea^ où l'oa a ootta aeulemaiit 
1^ tarmaa a&dés de i^ et dea puiasanoea supériaarea do &» 

z = p+pt+p'e+p'^e 

X = M + M'# + MV+ M'"f* 

i Y ;« N 4rN'* + NV + N'"«* 

Ces valeurs doivent être substituées dans les équations (i), 
aloa' de. i:éduire au plus petit aombce. les coefficifiiia indé* 
t€^i»miii4^9 et comme on. a 

ai on fait pour abréger 

t «ar mm'+nn'+^pp\ 

on «vira 

</• = r« + ait/ + &c. » VaF:R' + aK<+&a}. 

d'où résulte 




\ 

.( 3) 

JSt la substitution faite dans l'équation différentielle en ^ 
donnera 



,m i= **— — *r • m 



On obtiendra un senfiblable résultat des autres équations diffé- 
rentieltes.) ce qui donnera 

' = "('- ^ + 17) *"•'(' -67-) 
^-"G-^ff +!)*'«' ('-65=) 

IV. Au moyen de .ces coordoi>néea , on peut déterminer la 
longitude «t de la comète et sa latitude C, par les équations : 

De sorte qu'en «géuétaL, pour un temi^s qnelconc^e* qui n'ex- 
cédera pas certaines limites avant ou après l'époque choisie* 
on aura tes deux équations : 

» ( J T + -J-) . i*(i- r+--r) 



1 



(4) 
(3) / / «• ^ **'\ 

I mil r- + jTI 

y. Soil poar abréger : 

OT = M + M, m'=M'+/ 

et de plus y 

Si on mnltiplie les équations (3) par i H — ;j — —^^et qu'on 
néglige toujours les if" y elles deviendront ; 






(6) 



On peat simplifier encore ces équations en omettant le facteur 

1 + -^ qui multiplie ft', r , p'/ car la suite du calcul prouvera 

que l'omission peut être faite sans qu'il en résulte d'erreur dans 
les termes qu'on doit conserver. 




(5) 

VI. Supposons maintenant que les trois observations don- 
nées de la comète répondent successivement aux temps : 

tn sorte que l'époque * = o réponde à Tobservation moyenne. 

Soient 

les longitudes observées a*, a, a% 
et les latitudes V" ^ b y b\ 

Faisons de plus , pour abréger y 

tanga''—f''y tanga^fy tanga!=f (7) 

iangj^ ~ ^ ^^S ^ _ ^ng b' , - 

^ ^o — o > -^ — Si "~~ — r = g • 



cos a cas a cos a 



Si dans les équations (6) on fait successivement / = o , f =: — f 
t =: V, on aura les six équations suivante^ : 

(8) 

Vn. Après la substitution qui se fait immédiatement dea 
valeurs de » etdejt», ces six équations se réduiront à quatre 
contenant les quatre inconnues /» , /k', /,p', sous forme linéaire. 
n faut d'abord s'occuper de la résolution de ces quatre équa- 
tions. 

Dans le cas présent , où les équations (8) ne sont exitctes 



(6) 

qu'aux quantités près de Vùrùte t^, il fatulrti Mttettre > -dafna le 
calcul de l'élimination , tous fes tem^s qui fireroietft de cet 
ordre ou d'un ordre ftupérieur, attention qui contribuera a 
simplifier les résultats. On pourroit aussi établir les formes 
des diverses inconnues (a, fi'j f\ p, en les représentant chacune 
par A+ B9 + Cr+X)fl*+fi«r+ F«'»4*^c.;bn teconnoîtroit 
bientôt par la substitution^ que ^^s qhatrfe inoènmues ne con- 
tiennent aucun terme constant indépendant de 9 et S'; que 
la quantité ft en particulier, ne contient pas même les termes du 
premier degré y et qu'Ole doit étve ^visible par 94', idis sorte 
qu'on doit faire fc ±= A<«'+ B «•«'+ Cfl»\ On trouvera yareilie- 
ment que les trois autres inconnues «doivent être de la forme 
Ad + B4'+ G0*+ D 9 r+EB'% sans aller au-delà , afin dé ne 
pas introduire dans Jes équations ( 8 ) des termes du quatrième 
ordre; et cette forme justifie l'omission que nous avons faite 

du facteur ^+ (13 dans les équations (6)< On pourroit donc 

par les •coëficietis indéterminés , efibetoer la résolation des 
équation» (6); mais oa y parviendra auasi finilement par les 
moytns ordinaires , et voici le résultat ^ ce odcul présenté 
S<Mis la forme la plus simple dont il paroît susceptible. 

VIII. Soit pour abréger y 

Cette quhaltite qui se présente d'abord dans le résultat de l'éli- 
minitkii mérite iMe atteiltion particulière : elle peut se mettre 
sous la forme 

(/'-/) ig'-r)'-if-r) ig'-g.h 

et alors on voit aisémfetit qu^te est de rdrdre W (<+»'), 
c'est - à - dir^ du troisième ordre , et cette extrême petitesse 
empêche le plus souvent qu'elle ne soit déterminée assez exac* 
tement par les données de rol)servation. Soit de plus ^ 



(T) 

M (fg'-fY) + N (««-/r') = B 
M if^-rif) + N (^-^«^> ^ B' 

M' ifg—fg). + N* (^'-^) « I^ 

M'U'ir-fg) + N'(^-^) = E'j 
quantités entre lesquelles on a les relations 

B« + B + B'=:Ma E"» + E + E' =5 M'a. (lo) 

/•B* +/B +fB = Na /•£• +/E+/'E' =? N'A 
^B» + ^6 + 5^-6'= o ^E° + ^E+6r'E'= oj 

les valeurs réduites de kos q«atre inconnues feront : 

S^à. ©A * ' 

/ = î:ni?!±il''(OB»— a'B') + "^^~*> (tE^^|fyN 

m 

p=^ •-:=|^^VB'^«VB')+îl^(fl^iH._y^.j.). 

IX. Avec cea valeurs » on calculera, celles des six coëffidena 
m , ajp^niV'^VpV- ^ 4**^ se fera par les éq^oations 

at^f* + M m!^ i^+lilf 

»=«/i*+N »(5=»'+N' (i»> 

et ces six valeurs ne contieiidront que les deux inccmnoea 



(8) 

^et « sous forme linéaire. Calcalant ensuite la valeur de k par 
la formule 

k = mm + nn + pp' ^ (i3) 

et substituant cette valeur ainsi que celle de • = 3- ^ ^ 

r^ R. 

dans l'équation 

on aura une première équation entre les quantités ^et r, dans 
laquelle ^ ne montera qu'au second degré. 
Enfin on a de plus les deux équations 

« 

«'•+ n" + p'*= - — -, (16) 

r a ^ ' 

dans lesquelles il faudra faire de semblables substitutions pour 
obtenir les équations finales auxquelles se réduit la solution 
du problême. • 

L'équation (16) appartient à une orbite elliptique dont le 
demi-grand axe = a 3 mats comme on suppose ordinairement 

que l'orbite des comètes est parabolique , on pourra faire - = o , 

a - . 

ce qui donnera une équation de plus que d'inconnues. Et cette 
circonstance pourra offrir diverses combinaisons dont on pro« 
fitera pour rendre moins difficile la détermination des incon- 
nues r et ^., 

Nous donnerons ci-après le moyen de vaincre ces difficultés 
et de parvenir à des résultats dont la pratique puisse s'accom- 
moder. Mais avant tout^ nous examinerons , avec tout le détail 
nécessaire, le cas où l'on a 4^= 0; c'est-à-dire où les trois 
observations sont faites à des intervalles de temps égaux. Cette 
condition apporte dans les calculs une jgrande* simplification , 
et on peut toujours y satisfaire , en interpolant les observations 
faites à peia de distance les unes des autres* 



\ 



j 



(9) 

X. Soit donc 9'= 9 , et les formules (41) se réduiront aux 



soivantes: 

M = — B 

SA 

Or j'observe qu'on peut omettre les termes aflFectés de 0* daii5 
les valeurs de / «^ /.. En e£fet^ dans les équations (6) on a 

négligé , par rapport à ft', la partie -r-j y parce que le premier 

terme de fc' étant de la forme A , cette partie seroit devenue 

^-7 ; et multipliée par 9 ^ comme Test ^a dans les équations (6)1 

elle seroit montée au quatrième degré; ainsi elle a dû. être exclue 
des équations (6) qui ne contiennent pas les termes a£Eectés 
de t^ ou 0^. D'un autre côté cependant , en considérant comoif 
très-petit 9 le coë£Bcient A qui affecte doit être très-grand , 
afin que A0 ait la Valeur finie qui convient a la quantité ft'; 
de sorte que A9^ désigne réellement une quantité de l'ordre 9\ 
Ainsi en conservant dans les formules précédentes les termes 
affectés de 0% qui ne peuvent être beaucoup augmentés par 
leurs coëfficiens , on n'ajoute pas à Texactilnde de ces formules 
dans lesquelles des quantités du même ordre ont été omises. 
Donc il faut supprimer entièrement ces termes , et on aura lès 
valeurs très-simples : 

f^'-^CB-'-B') . (17) 



• • • 



( 19 ) 
/«^(/«B'-fB') 

SA 

dans lesquelles rinconnae C ^^ ^^ retroavc plus. 

XL Avant d'aller plus loin, il convient de âéterminer les 
quantités M, N, M', N^, qui dépendent du mouvement de la 
terre. Et d'abord y si on appelle Â la longitude de la terre, ou 
celle du soleil augmentée de 180% au moment de la seconde 
observation , on aura 

M = R cas A , N =s R « m A. 
Appelons enpore pour un temps quelconque / : 

f Texcentricité de Torbile terrestre^ 
* la longitude héliocen trique de la terre t 
y la longitude de l'apbélié , 
V ou R le rayon vecteur , 

on aura par les propriétés du mouvement elliptique et eh ^ 
regardant > comme constant , 



1 — •€»«(♦ — y) 



R j. u\ 

dt " R' ' dt l/(i — •') * 

Orona X = Rcotf«^ y=Raî/i«; de-là on déduira les 

valeurs de -~- et -^y lesquelles en faisant f=o ou * ==A^ 

dt dt 

deviendront celles des coefficiens M' et N' : on aura donc 
, — sin A + f sin y „, cosA — %cosy 

XII* Si au lieu de la constante y on introduit l'anomalie 



V 



(11 ) 

vraie **• s» A •— > , ce qui donne > == A — t^ j on aura , en 
développant ces expressions jusqu'aux §* inclusivement, 

N' == co* A ( i — -« « C05 **• + 7 é» ) — • ain Asin ^. 
Par un semblable développement on a 

« 

i 



Th sorte qu'on peut écrire ainsi les valeurs de M' et N' ^ 
M' = — sinA f — 5^— J ! — • Vi/ï ^ C05 A 

d'où résulte i 



(«9) 



(ao) 



(91) 



MM'+NN'= — ««n^.R. 

On- conserve dans ces expressions le quarré t^ dé Texcentricité', 
quoique^ ce terme soit très -petit , parce qu'il ne complique 
presque pas le calcul, et que si on l'otnettoit entièrement , la 
valeur de R se réduiroit a i + t cos **•, et ne correspon droit 
pins assez exactement au logarithme donné par les éphémé* 
rides. On peut supposer pour les temps voisins de celui-ci, 
f = 0,01679, alors on a •*= o.oooaSig, et 

^^\ '"r'' / ~^^f^%'^ o.Qooo6ia, 

A l'égard du terme t sin Y, on le calculera en prenant dans 
les tables le lieu de Ifaphélie ^ qqi CîSt très-peu variable, d'où 
l'on conclura l'anomalie f = A — >. On peut encore déter- 
miner la^quantité tsin^ avec une précision suffîsantp par la 
seule' valeur connjie du logarithme de R; car tirant de ce 



( i» 

logarithme le nombre -^ s= i + « , on aura » co* ^ = «• — tt , 
etde-îà 

• «»* ï==fc !/[(. + «•— I,) (, — ,» + „)]. (92) 

Quant aa signe de celte quantité, il sera positif depuis le pre- 
mier juillet jusqu'au premier janvier ^ et négatif dans les six 
autres mois. 



i 



* «t 



1 - - •— f 
Enfin on peut encore observer que les quantités — ^ — et 

êsin'ir^ qui entreront comme coëflSiciens dans nos formules, 
sont liées entr'elles par cette relation : 



(^) 



9 a 



Xni. Il s'agit maintenant de former les coefficiens des équa- 
tions à résoudre, de la manière la plus simple et la plus com- 
mode pour la pratique, PouY cela, il convient de remettre à la 
plaôe de /, g^f^y g^y &c. , leurs valeurs en a , b^a^yb^y &c. 
données immédiates des observations» Cette substitution étant 
fidte d'abord dans la quantité A , on prendra pour abréger, 

D = tangb' sin{a — a*) + iang b^ sin {a ^ a) 

+ tangbsin(a^ — a'), (a3) 

~D 

et on aura a = -,. Faisant ensuite les mêmes 

cas ar cas a cos a 

« 

substitutions et celles de M = R cosK , N = VisinK , dans les 
valeurs de B% B , B', on aura 

— = — =r — Itang b sm (A -^ a*) — tang 6* sm (A — a)] 

A XI 

IT IXcosq 

A ~ D 



[tang b^ sin (A — a) — tang b sin (A — a*)] (aij 



B* Rco5a*_. . . •A • /\ .# • /A \-> 

— = — =r — itangb 9in (A — a ) — tangb sm (A— a }J, 

A Jjr 






(i5) 

Soit encore pour abréger y 

Ç = tangb' sin{A—a') — tangb'sin{k—a')j (a5) 

la première des éqafttions (17) donnera ft2= —z^Ccoâa; or 

par Tes triangles rectangles TdC, TKG (fig. i), on a TK = p 00s i, 
et TG ou fi = p Cas b cas a } donc en introduisant p à la place 
de M, on aura 

_ Rijg'C _ Rfl'C (}:^_2\ f ex 

' ~ aD C05 6 "^ aD cas b Vr' R'A ^^^ 

Substituant de même l'expression de fi en p , dans les valeurs 
de i7i| n^p^ des équations (ia}/on aura 

m = p coJ^ a cas 6 + R ^^^^ A 

n zsz f sin a com 6 -4* R êin A ^^7) 

p S= p 3m 6 y ' . 

^l réquati&n in*+ »*+p'= r* deviendra 

r»=R*+ aRpeej ftcatf (A — a) + p*. 
Cette éqnation seioit donnée immédiatement par le triangle 
CST; car en appelant c Fangle entre le soleil et la comète 1 
les côtés compris sont R et p , et le côté opposé r, ce qui donne 

r^ = R* — aR p CM c + pV 
Or Pangtê c a pour mesure l'hypoténuse d'un triangle sphé* 
rique rectangle dont les côtés sont b et 180^— *A+a, on A 
donc cos c = — coB b cos (A — a). 

XIV. Les coëffictensC et D se déduisent immédiatement des 
données du problème par les formules (2i3) et (95) ; si ensuite 
on fait 

cos c = — cos b cos ( A — a) , A = —p: r, (a8) 

on aura les deux équations 



t 



( 14 .) 

dans leaquelles il n'y a d'incoimaes- que les deux côtés' r et p 
du triangle SCT. 

Si on élioBine p de ces deux équations , on aura une équation 
en r du huitième degré , mais qui sera divisible par r^^^R , et 
s'abaissera immédiatement aa*septième. Cest le résultat auquel 
sont parvenus tous ceux qui, ont traité avec succès le problâme 
des comètes ; et ce résultat est d'autant plus remarquabjle , que 
les équations (39 ) se déduisent assez facUement des données 
de l'observation , et qu'elles ont lieu sans supposer que l'orbite 
de la comète soit une parabole. . . • 

XV. Comme les quantités r et p sont essentiellement poisi- 
tives, il n'y aura d'admissibles, parmi les diverses solutions 
deef équations ( 29 ) , que celles qui donneront r et p positives. 
Or en construisant les deux courbes de genre hyperbolique 
qui représentent ces équations, et examinant les diverses iuievr 
sections dont elles siont susceptibles , on parvient à cette con- 
clusion générale sur le ifombre de solutions utiles que ces 
équations comportent:. ;'. . 

, Si .3 h' Qos c est positif et plus grand que R^, les équation, (99) 
admettront une solution et n'en admettront, qu' une ; dans les 
autres cas , ces équations auront deux solutions, ou n'in auront 
aucune. • 

XVI. A l'inspection de l'équation P^^V^J*"^)? on voit 
que le signe h doit être le même que celui de 3- — —^. Donc 

r R 

si h est positive, on aura r<R, c'est-à-dire que la comète 
sera moins éloign^^e du soleil que n'est la terre. Au contraire, 
si h est négatif, on aura r>R. Dans le cas où l'on aurbit 
D = o , A seroit infini et on auroit r7= R, c'est à dire que la 
comète et la terre seroient à égale distance dû soleil. Ce cas est , 
d'après la remarque de Lambert , celui où les trois lieux 




\ 



( t5 ) 

lipparens de la comète seroietit situés dans le plan d'un même 
grand cercle; car pour que trois points dont les loûgitudes sotif 
a% a , a', et Jea latitudes b%' b , b\ respecliveoieal, soient situés 
dans le plan d'un même grand cercle ^ il £aut qu'on ait Téqualioa 

o = tàng V s in (a — a**) '+ tang &* sin (a' — a) 4- tang b sin (a^-^a). 

Donc lorsqu'on aura D = o, ou seulement D très-petit, ow 

pourra en conclure ^ ou exactement, ou au moins par approxi^ 

mation^ 

r = R et p = aR ces c. 

» 

XVII. C'est une circonstance peu favorable que celle où l'on 
tombe sur des valeurs de D trèsrpetites; car alors une erreur 
très-possible d'une ou de deux minutes sur quelqu'une de$ 
quantités données par l'observation, pourroit changer dans une 
proportion notable la valeur de D, ou même lui donner un 
signe contraire à celui qu'elle doit avoir. Pans ce cas , il n'y a 
aucune conclusion certaine à tirer des équations (t»9), sinoii 
qu'on a à-peu-près*r = R et f==uKcosc; on se dispensera 
donc alors de chercher les jracines de ces équations , et oh aura 
recours .a>ix autres équations qui seront. données ci-après. 

XVIII. Dans les autres cas , la résolution des équations (29) 
s'effectuera assez facilement par les fausses positions j mais pour 
cela il est bon de connoître d'avance les limites de r et de/. 

Lorsque h est positif on a r < R j mais la seconde des équa- 
tions (ag) donne r'= R*5m* c + (p— Rco5 c)*^} donc on A 
en même temps r>Rw/i c. Ces deux, limites auxquelles .qp 
» peut joindre celles de p, savoir p>o et p<all coi c , sont 
suffisantes pour diriger le calcul et rendre la résolution assez 
prompte. 

Si h est négatif et cosc positif , soit A =: — -i; ayant déjà 
r > R , on doit avoir par conséquent p > â R cos c; mais de 

l'équation -^ ~ vT^ — - y ^^ ^^^ f < fTh h donc on aura 



( i6 ) 

r<v/ (R*— — ^ ''pir ■'^ **"^* observer que si oti 

n'a voit pas î > 2 R^ cos Cy les deux limites de p seroient incom- 
patibles, et il n'y aaroit pas dç sol a tien. 
Enfin si h est négatif, ainsi que cos e; soit /Is» — î, et 



cas e^= — /, on aura toujours p <! 51 , et l'autre limite sera 

MX 

simplement p>o : celles de r seront en même temps r>R 

^ y/«. 2 Ri/ i*\ • 

etr<v/(^R-+..^^ + -,| 

XIX. Il faut maintenant faire la substitution des valeurs 
de m'y n'y p' y dans l'équation (16) : or en vertu des for- 
mules (13), (17) et (a^), on a 

^ , R« • f cos a*" ton^ h sin ( A-— a') — co5 a* tang V sin (A — a)| 

3D .( +C05 a' Umg b sin (A— a*) — cos al tang V sin (A — a)) 

. . • . (3o) 

^ _^ Rwdf sin a!" tang b sin (A ^-^ a') ^-^ sin a"* tang b' sin {A — a)) 

3D l+ima'ftiiigrd^J's (A— a*) — 5wia'la;^i*«w(A — a)j 

i 

, R»9.f tangb*tangbsin(A'^a)'^tangb''tangb'sin(A'^a)] 

^ aD [ + te/î^&' tangb sin{^A — a*) — tang b' tang b"". sin {A — a) j 

Appelons F, G, H les quantités qui multiplient — ^r- dans les 

Vadeursde m', n\p' respectivement ; nous aurons en substituant 
aussi lès valeurs de M' et N' : 

R«« 
m= — «7» A I — =r- — I — • sin '^cosA + — =r- F 



r = — «7» A f — g — j — • su 



ri =s cos 



aD 
R<id 



• « 



A f "—^ — J — tsm^ sinA'\ jÇ^j 



d'où résulte 



/" 



( 1? ). 



Donc m. pbar abréger «ncbre les «xpressioiul on &it 



f 



(3i) 
Q=: fang b tin (t^^c^^^"^ ~A»»ir**«>ç(A^a»)«m(A7-a) 

— teiy i»CM (A— a' ) CIA (A— a) 



H=s — a «ai,^ je HB^i* rt» (A— a) 4- #«1^6 toiv i» ,i„ (A_^) 

+ <âiy6<(M^i' MI»(A — o^). 

Le» troia ^luaatUé» w', a', p\ poorrant «toe auset «otis oèttè 
Sua»: . ,,...■.. ..•••, 



et leur Bubctitution dans i'éqnatlon (i 6) donnera 



la,»'/, — i,-^ , X 



XX. 



3 



J 



( >8 ) 
- =: o, et réqoation précédeble derra être oombiaée a,T4^ 
réqnatioa 

d'où, l'on Yoit qae la rés vivante- qa r ne serfi qae du sixième 
degré. 
Cette éqoatipn aura toujours Vue raciae réelle potitiTe ymti 

en &isant »s=x, - =^^, si. on constroit' les deux par^boks; 

qui o&t pDor équations 

■(.■■■■ ■- ■"■ ■-■■". ■■ .■ 

La première sera une parabole ordinaire , ayant son axe 
dirigé jdans le sens d«s drdoanéqs. Cette «cnirbé est aitnée tont^ 
entière du côté des y posîtiTes, et ne coupe point la ligne de» 
abscisses , parce que le second membre de TéquAtiorr ( 33) formé 
qe la somme de trois quarrés , ne sauroit,deyenir nul et a 
ses facteurs imaginùres. D'ailleurs lorsque x := o , l'ordonnée 

i;"'^-^T. opà-pen-prè^ j-s^f ' y 

1^X41 seconde courbe est une parabole cubique don( nue branche 
à'étend à l'infini dans l'angle des s et y positives, et l'autre' 
dans l'angle des x ety négatives. Lorsque « ^ , l'ordonn^ 

y^-^ft où à-pen-près j{^== i; ainsi dans cette ^ourtj^, les 

ordonilées sont plus grandeaVers l'origine des abàcissOs que 
celles de la première parabole^ Maw comme à une grande dis- 
tance de l'origine c^esî le eontraire qui a lieu , il s'ensuit qu'il 
y aura nécessairement une intersection ; et il ne peut y en avoir 
qa'wu»^ 1/vr&) que la jiarabcde' ordinaire oppoH w'coiivexité à 



(19) 

:la]ijgiia:dètidbMiaiés , Ifthdis que la^Mlkde cbbiqaô-IuitoppeM 

•«a. cmictTilé ( f^M fig. 4 ). 

- Le syatAme ides équations (S3) et (34) estidoncpréférablê^ 

uceloi des: équations ( ag ) , non-^seolement parce que le degré de 

l'équation ^ôale est moindre, .mais encore parcje que la soio- 

4ion^de celles<-là est toujpurs unique | lèt ^se préseiite aucune 

jambiguité. 

n fioLut observer néanmoins que l'équation (33) est eu)et(e 
an même inconvénient que la première des équations (sg)» 
tdaos ie cas où D est très-petit , parce qu'alors tes èrreu rades 
observations ont une trop ^ande influence sur la valeur de X^ 
^qu'dixisii la résolution des équations (33) et {^34) ^oftrroit 
mmduire à des résultats défectueux. '! * 

* • • » I 

. XXI. Pour obvier à cet inconvénient » il fant éviter tonl 
•emploi du coefficient Dy et on en a lieureoseoient .la .fiu)iUté«| 
' puisque le problème des comètes y dans la supposition d'une 
vorbite parabolique y o£fre une équation dç plus ^ue dHnoonilues. 

* Par Téquation (a6), on a — =r* = — rry^^ et celte valeur 

étant substituée dans l'équation ( 33 ) > on aura^ en suivant: toft^ 

•f 

|ours l'hypothèse parabolique y c'est-à-jire eu faisant -r^Of 

• •• 

a a 

- '^ ^ , '.* ' ~ (35) 

p* cas b . \ * 

équation qui ilevra être combinée «Tec Téquatioii 

;'=R•-rp(aR«a«0 + ^^ (3^) 

afin d'en tireur les valeara de r et f • 

Si ou éliminoit r ou f 4^ ces deux équaliona y l'équation 
finale, nie monteront encore qu'au sixième degré. D'ailleurs les 
cpc^cieiis de^cea éqnatiank «c; diédoÎMOir ^▼m ^^^ exactitude 






) 



-tfiiSMitttf, ûeê données do l'obtenratiotL I^ ooèfficiéntC qiii 
sert de diviseur ^ est celui qu'on coQiioitraaTOO le fins. «le fié^ 
ciéion $ ctkT en regardant 9 comme une quantité éa premier 
4irdre 9 C^»l du pnmier ordre aussi ^el les Iroîsaalpèa P9Q9 H 
-do seconde II parent donc qàe la c<raibsaaâ86B>des éîtuatimu (35^ 
#t (36) réntlit tous ica aTUitages qà7o« peut désirer dana là 
solution du problème des comètes ^ et c'est à celte combinaison 
-qn'iJ connrient de s'arrittkr définitiTumetit ... 

XXIL On peut prouver encore que les équations (S5) et (3Q 
aiMnl toujouta nne aololion en noMbrea poaâlifsà. 
^ Su cffel y ai on décrit les denic xHKtrbca qdi tépohâent 'h\tsn 
équations, en prenant p pxmr Tabseisse et r pour Toidoiinéei 
le lien de l'équation (36) sera une hyperbole équilatère dont le 
lùmUte est sur la ligue des ubseisBes à ia dhtaHM A oé^cV^ et 
idunt ludetnfaaicie, dirigé euitutit les eidemiéts «a R sét o{r&yek 

» Uéfjfitâlon (34)^ qu'on pourrait réduira k la forriie 
y(af^'f'a^):=:a^hy appartient à une courbe qui ^ comme la 
concboïde supérieure , a pour asymptote la ligne des abscisses , 
^ aasl située toute entière d'un même côté de cette ligué , sâf 
plus grande, ordonnée étant h. iLorsque p=so., rovdonnée 
de lliyperbole équilatère =R , Celle de la seconde courbe 

= — ^naJB..à peu«^près; dàxic'paisi[]nb la première s'éloigne 

à l'Infini de l'axe, tandis que l'autre a'en rapproche de plus en 
plus à mesure que 4'abscisBe est plus grande, il s'ensuit que les 
deux courbes auront nécessairement une intersection dans le 
sens des abscisses positives et n'en auront qu'une. Ainsi la 
résolution des éqrfatiotis (3f?) et (3fey*auVa encore l'avantage 
d'être toujours possible et die Yi'effrir iMMile ariibtguité.' -> 
' i\>ur avoir les Utnites de r, il fiMitofe^hver que les-secends 
membres de» équations (^^ ^ (96) ne f>éuv«nt s^ rédati^eu 
zéro 9 et qu'ils ont tooë 4eux un minitmèth^ Le Minimum de ^ , 



i 



((tli) 

m«iidN3« à» réqiMtioil (Je), et il est clair (^a'oq. «ma i Ufois 

r>R^iii c et r< -=-. Au moyeu de ces limites, il sera tod- 

jours assez facile de résou^lre les éqnatioas dont il s'agif par les 
l&usses ^K^iiiôns. 






XXÏII. L'applitoUon de nos formules exige ^ue l^t^eoips 0^» 
donqé en jours et fractions de iour , temps moyen ^ soit con?e;5ti 
en arc du moyen mouTement du soleil ^ cet arc étant lui-même 
rédiiit en parties du rayon t. Fôhi faire cette double réduction, 

. . :.*.' • l '. fit; m'é*d- -<« ' JiN.J'/ '••;IM * * - Il > ) 



I . 
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dra ajouter au logarithme de , compté en jours du *tèttf^ 
moyen, le logarithme con)|tant.*3|33$^SK^u 

En faisant abstraction des erreurs des observations^ les coëffi*» 
cWs ^i'^TëéX'^êtéÊdVê^ ëékiHê 

qu'il faut choisir trois observations faites dans un intëlr^iftSIè 
de^mps peu. considérable , et^qoj n'excède pas quinze à vingt 

joura. MW3 il, IIP fj^^l- pas noi,^ 4«*«,aflf jç?^o^?.^»aijpi|ft ^^^ 

trop T^jB^i^h^s ,. j^i^pçjqu'a^or^ ^p m^^^^S «I^^^M]»ft? 
Jel^iÇP^ète,if«^pi^txop,petit^ çt Jçs çr^eft/ft^Ef P>?ÇfiW^fiW 
auxaient,|înc^ UQJf^m^^^^^^ ^^.^9f?^^^ n;;iî'.r 




qn^il y aura toujours liei^ à les porrig^er par le moyen 
vation»^^ plM éloî^fiées iBir^l<Q» , W pp4irva/è||P ^^MUtiw-Jc; pkA 
souvent à chercher une valeur approchée de la distance péri* 
Blêlîe , et de l'instant j^u pqssajge, atf ^erih)llie^,çiur aveerl^ con- 
noissance approchée'^de ces deux élénlené ^ on est ep état de 



î 



■'pvouéder «BVftkafi néoestaive'i^oôrr CtbUtair uiiè «ni^ite )ebirigée. 
'Or sl.oii4ippelle-n la ^tan6tt pénhâie, -et.qu'oè oatcaleX 
par la foripale 

comme cette quantité est égale à —7— , si elle est n^gaUtvii^i 

on saura qoe la comète s'approche de son périhélie ; et si elle' 
est positive y on saura qu'elle s'en éloigne, pe plus , au pioyen 
Ûei^ QXi, aura la dbtance périhélie 

ft=.rr-i^. (37J 

Cette distance étant connue, on déterminera l'anomalie vraie 4 
de là ronlëte au mopMot de^la seoofide ' observatioii j par h 
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jEnfiii l^fmoiQaUe étant connue t:Pi|chcancheaa dam la tabl^ 
4ifi^ cpni^lfis )p temps T qui répond fi çottp «nomiUe ; liiis^nt 

ifHlsaitfr :. , 

• * = niT, (39) 

'bu àuf a le temps i employé paif'Ia comète k parcourir Tano» 
li^alië' ^«^1 Çè^* temps étant ajouté à Vépoqué de l'o^bservâtiofi 
inbyeiine, 'ai ta cèniéte avance vers son périhélie, ou en étaiit 
retranchéV^^ ^Uo s'en éloigne , dbdnera Pinstant du passage au 
péri)iéUe. 
,'l!V)ut se réduit ïonp à déterminer la quantité t . dont la 

§*a^eur et lé signe sont .également nécessaires &, considérer, 
r'en stibstittkant les.Taleiirs données par les équations (3â)« 
on trouve 

^ «s f cothêh (A«-«) \J-7ô"^ "~ ^ B Kp ***** 1^ ** 
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i XXV. Maiâ «i ob vent ^étonmner. à la fok tous les élémeas 
âe l'orbite , il sera inutile de chercher partseUemeat la valeuf ^ 
de t^ et il faudra calculf^r les Tueurs des. six quantités tn^ n^ 
p , ïn\ n\ p\ par les équations ( aj) et (32). 

Ces Tsleurs étant connues , on çn déduira itnmédiatei^enl;^ 
celles des trois quantités mu'— m'i», mp'^—m'p^ np''^n'pn^ 

La wemière étant égale à -^r^^ « ^ . fera «omMltre mar ^n 

siffue si le mouvement de la comète est dirept ou s il est rétcor; 
grade. JQans le. premier cas, mn'-^ m'n doit être positif^ et 
dans le second il sera négati£ 

Appelons maintenant I ^inclinaison de l^rbite , ^t S la Icp- 
gltude de celui des deux nœuds que la comète deyfince ^[ans 
l'ordre des signes. On aura , par les propriétés du mouyéinêné 
pars)K>lique , les trois équations : , 

m»'^— » m'n = db cas I j/în 
i I mp'.i-^t^pifp =^ i±:sinIcMBf/ûn' ^ ' (4») 

np' — n'p = ± «m I sin S ^/sTn; ! : . 1 .J ^ j 

pam cea {iMmul^^^qni doÂy^nt; V^ijoi^rf .^ounçiç.^^ et 
l/iril positif, il ^aut prendre les pgnes ambij[us du second 
membre tous ^^dsitiVeiAènt ; lorsque >7lh^—^mn est positif ^'iX 
tous négativement ^ lorsque m »"— m^n est négatif. ' ^ '/ 
' Là longitude lïén pàrifcdlier sera dëterniinée par la idtiniAe 

M»i# eett9 Icuteule, qui convient égedetomift mm BMfjibsc:Bi>€l 
i8o^+ S y donne la position de la ligne des nopnds^.et n'indique 
pas Fun dés nœuds plutôt que Pàulre/Oii'àcKèVerà'de détér^ 
miner S en calculant l'inclinaison I par la formule l 

>et prenant des deux jfipxm de;^cejlpqi;i /»^i/4W«L^^<Vrm 
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et on lèveroît toujours dé ta même lÀahière riadlétermination 
de S. ' . 

X^ iMfM^f^timiî é^t^iàinj^ par la vàhtir dé 8 j >anti iW imœmd 
ascendant j^ si le moupfmeniest dîrecf et ta latitude boreçle, ou 
si lé mouvement est rétrograda et ta latitude austrçïe. Uàns les 
aeiùc autres cas, ce sera le nœiid descendant ^ et il faudra 
ajoute^ 180^ à S pour at^oir le' lieu du nœud ascendant» 

t^8 àagleè S et I étant détermines , ,on aurait la d^^tance 
penhéhe It par la formule 



■'•^ ... 



/ 



^^(mnV^r^- ' = ' " : ^ tin 

0m:]^urroit anaâ ]|i rtmoTter indépendaMm^ii de ces angles 
par la formuler ^ . ; .; '. 1 .... .^ v. — ' . 

^p jjj^jiltQ ^xw^mnieijt ^eff é^if alion» (4i) j et il est aké de 
voir qiie ce^tç^e^i^^ valeap s'accprde avec la formule (37). 

décompose en ces quatre quq,rrés : 

Doâaifc.aoïaBiw dcftrttbiiidanûcia qèafrés ^ (mP .4.'n(*-f /^) 

XXVI. Avec la dist{ince périhélie n et le rayon vecteur r, 
^.^Iculera Paupaiabè vi?aèe4^,de4a«>]ii^.au moment de la 
conde observation^ et ide^TàTinalant du passage au périhélie, 
^ dtfoki\i\BipHi|ué dansli«rt;^XÎX3fV^ ' ' ^ '- '^ 



Enfin povr^ avoir le lien du périhélie ^ on calculera la lon« 
gitad'e li^ïrôcetitriqae P Ao là comité par là é}jrîàidé ''^ '^ 






« r f 

tang9 = -^y (*9 

ayant soin de choisir entité' les deux angles V et ip4- iBo^ qui 
conviennerit également à cette tangente , celui qui rend sinf de 
même signe que Oy on) cas ^ de même signe qde m. ' 

9 étant ainsi connu, on trouvera toujours <|ué^-^S éA 
positif et plus petit qbe iBo^y conformément a la supposition 
£ûte que S est Ja longitude de celui dès deux' noeuds qui est 
moins avancé que la comète dans..Vprdre des signes. Cela poséi 
Hcttt ^ la distà]ii3è''dé'là ééibètè aà ^qsttl'dont la 4<Migildde èk 8\ 
on calculera ^ par la formule 



■ •.;. ' j; : . : ; ' ^.. . ^ r." . 



et le lien dupérihétie 9ur Torbitcseta. 

. <*' \ . i \ - ' ■ • > .î _ i » - « f ■ > V * ' ; i 

S + ^ + 4 • d lacoinite mwcfae iwrs son périliéli|e y \ 
04 S + ^—4 si elle y a déîÀP(S8^ 

XXVIL Pouré viter la peipe de cherchéir aiÛeurs la démons- 
tration des formules (4i), nous croyônTdevôir rajouter ici 

Des équatioris cKfférentielles (i) on aeduit immédiatement 
les iQté^aleé sqivànfteB^ qài. ésparpneni.qiif) Ibfr'airèflr décHteâ 
dans les diverses projections de l'orbite sont pnDporttckimdUes 
au temps * ■ '.-:''. '» 

' xdy^r-ydx =<Xdt 
xdz-^.zdxi =; Cdâ (48) 

on en tifis d'abord Téquation ..... 



féiMnl 
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des noeod» - :^ te/i^ 3 , nmi . P»! *!ïr* C'' =; C" <«»« ^ D» 

rayon i décrivons dans le pian des « et^ Tare indéfini « dont 
l'origine soit sar l'axe des.a7;,IeÀ coordonnées du plan qai 
répondent à ^'extrémité jlq c^tarc seront; ;,,.[., ^^ . . 

X = cos s ^ y ^ sm > • ^ ^ = 7^7 5"» * '^-" tït (^os s. 

Ou attriT done ra«si ••' •••!■ i-' •' i* • .:■ .''•••> ; ■"• ■ 

,*t à^wM l^emph ^^^^^rg^ 4«ivi«it.M *img»^ de l' 

son ou to>ijr I : on a donc 

donc C = C co« S <a»^T et C" ^é sin HiàÀtl':- 

CeU^otfè^, kè ^qtiàlfôli» |48 )'rïi>pdt»t«6S «^I^ot^éx)& l'on a 
»=i: m, ^ = », &c. deVjenient-' * ::;;•. ; 

. r ,'mjp — fft'p = c cosStangk 

et iVnê reBté*à dëterfnikier q«i9:£'.iQr dam r^qoAkiahi ^16) icpii 
çétLt«^rii«^4iiqsIr . :.**!• '' • ^ <n. ij./- i ; • .. .-.v , . ' - % 

le premier membre «i > ■ - ^jm ^ ;y*\ multipliant de part et 

d'autre par r% et mettant à laplace de -^ sa valeur ï, oit èfi 

déduit 

~ _ ^^^^ ^^ ^ * 



» ' 

▼ 
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( ?7) 

(jf^ a^ cas I C04S 1 ' 

coi*^ I \i a tt. — — L ïjfôric on a eh' 'générât , ' ' 

et ei>|iarti)ià1iei^V'for^iie 1\>c%it6 est ptnbbliqte, 

ce qui donne les formuler /4i). 

' JC^tli^l^a^tilaKôn^^ûe Mous venons ^^S^Velo^rpwif 
fë càÀ'bliléï'bbsrérValfôfi^^^^ Ar^dé* tnJèrv^les kle tèups 

ègau:<; paiÎQih drrèlr îèrtife TexftDtif^ M la stài|»Iîibîté i|«'oa 
peut désirer dans une question dont la difficulté est génénd^ 
meht^tedbuniié.'Si par dès liifèdiMtancës ^\i ftvohibles cette 
solutidii hë conduit pas à deé rësnltirt^ëiiMsam'mènt' exacts , il 
faudra kietï' prendre ii! Fimperfectibil: de^ <^bser?atio«is et non 
àaxil<égltgèÀces quéiid» ri<ms s<$tntf)e9 permises* dans la réduc^ 
tion des fôrtnnle^. ~fin i^aîli vofadrbit-on pottsaer' plus loits 
féxâétitbfi^ de des Ibtmèles i^-ik i^di ne pbftrroit ^ fiiire qœ 
Mi: ïlfes èktcôls'' tro^ ^oc>mpti(p}ésr-piFnr-<M^e de quelqn'uftilité 
dans la pratique ) ; les erreurs des^ èfbsèmttions domin^tfoiènt 
iéjttjiburs et empêcheroient lësHrésyâtals de passer un certain 
^égré^^>9nf xîuati9m D'i^iHeurSi on .wit «ssçz que trois obser^ 
«atKi^ÉHpiiàcip »p€f^ de di«lfACç le^ iûlÇ& dfa.tjUtireSy ne peuvent 





< 58 ) 

pas détérbiider aïisëé eitaMe^eitt llDtÙfe-dVtt^^émMe/ Ut pfiM 
faiit topjotirs 'i^bthrir i des obserViBLtidns ^^8 é^igÊtééà^paaé 
r^cniSérles^ étëmens déterminer pér^^ 

tion. Il est donc iqulile de pousser le scrupule trop loin sur oelto 
première approximation^ et Pobiet doit être.' suflfisamment 
rempli par la méthode que nous^aVbns exposée. 

Lorsque les troispbservations données ne seront pas égaléttiënt 
distantes entr'ellés,. il>faÀadi»i,ca3|iiq]^ n^rt^^ravons dé)à dit, 
calculer par l'interpolation au nDuveau lieu qui soit ayec^eux 
des autres dahsia conditioii' ri^iiisi de Pégalité^^ldâitaocei 
Pans ce calcul d'interpc^lationj il ne faudra fiâire enlxer^que lea 
trois observations données / ou en général les toôis observa- 
tions connues les lpltis_propIiW da mom/ent pmir lequel on veut 
déterminer le noat*eau liem Tels sontles préceptes au moyen 
desquels on phftiettddfa toujpiv's nnesçlnJim.sjaffifi^vun ^pp^f P; 
cbée dans la pratique, . • _ ^ . . 

XXIX. Mais à considérer leprol>iéme sous un point ^ë*vue 






porqment^ ifnaly tiqiae :^ aqi)s nVM^iriQnf.^^aqpiU.qia'ini^ 
aient: oMre bfajeti si nc^ ^e liiiaions yaic |wr qppls inoy^ns 
on peut ]ob tenir une solution* générale indépeudanta de l^jfnp^ 
pabitioâ:de -0''.s7^0.. "...',: ./^ 

Dans 0etta yue , nous a\loa3 tepirçndr^ Ips jéqoatfons (\i} > ejÇ 
bo^is. observerons .d'a[bQ)r<l.*qaÇi les piêmes . raisons pf c jte^ijttelles 
lions Aypns supprimé Jes tei^qies i^Qecléade fl* <iUns je^ {prpauïea 
de. Tart. X^ doiv/eipit égalemeq^ faire omet^tre, dans les valeurs 
d)OM Qt ty l«g termes. c^eetés 4a fjic^eaç 6'y- ii'-fr j9'*f.i[I.Qe reste 
akjTs^ 4U|i« lQ$<àluatioq& (i.i) que ^eux sof^s d,e ter^fa.^ ïe»^ 
on» .^p(éf d^Jaclevr.,•-^^(fl'^fl)„le»,aJlt;PÇ8 affeçt^ .^^ 

le second facteur »(«*—») pourra iSti-e re^iréfleWférpif K^-i^f;)^ 
car )a différence de ces deux quAhtité«>/4éi «Àf^l'-M^yM 



\ 



(«9) 

ttonve de Fordie des quantités qu'on n droit de négliger. 
Il soit de-là qu'on peut mettre les équations (ii) sous cette 
forme: 



i 

aA 



/,'= •::V(«*'B--«yB')+-^(«iS'tE«-*V'B7 



Soit V xsz 3/+ y et t =s «' — y, c«6 éqintiom defiendront : 



M = 



M = 



y'» 



/>' = 



aA 
aA 



(B+ ^y^) 



•i» . 



* * 



[B'— B'+ !/(«•— E')]+'-Çyïiii- 
[B» + B' + fy(E« + E')] 



j; /•B«-/'B'+4/(/* £•-/'£')] + Ijr'N 

[y^B«H./:B'+ f y (/'E'+^E')] 



, [«^B'-«-'B^+fy(5^EWar^E')] 



fi) '.o ..''. j) !'■ ^ • -1 q ' '• •' ") '' ' ••' '■■' ' ';• * '■••'-• • 

dè$.!rAMi»ife-Bn+B:, /^B'.+/BS ««.4«iiiéB»j^Bl0»éqiia- 



«'_-* 



*■* 



"= \'£:'cB+-;/«r.:; 



r • » 

. î. • » 



» / 



4 

^'= ^[B'-B'+|y(E'-E')J 

a A 

I - ' ~* 7 ■ • ■ . ■ 



• '4 



» - f 



P 



. t/-B'-/'B>f/(/'E'-/'E')} 

SA 
2 A 

2 A 



+'^^iCB>TyB> 






. XXX. Maintenant la substitation dàs VéUursr àt f^ fj /'t'^c* 
en a^ b^ a\ &c. , donne, aairànt les dénominations . eipployécs 
dans les art. XIIi^^:3(];;X» } ! 



. * T 



B _ C 

— = H. t»«.a."sr 

A • - S 

B*— F „ F 
fB'^fV _ G 



s I 



" t 



^W-g'B' _ H 

A l'égard des antres (Quantités composées en £ , commor celles-ci 

le sont en B, [l lila1^bba«^iHBÉ* 4«>e'ï:4i^^^^ 

dans leqnel on met M' et N' à la place de M et N. Or en 

Btem >{ià]i0S( ^tttt MU» Itt twiites; a&âlôs^e; Si<fiit^«»iiAiplléi^ 
par la quantité très-petite /), on a M'stt:L^^'«ii|xéhrN^tt±cM«A, 



:(îi:) 

on encorçy par la hiènie raisox^ , M^ >= — B. ^m A j N' s^ R cos A. 
ObservoDS déplus , que cliaoune des quaiikitéa G , F, G, H » 
considérée comme fonction de A, f«t de la forint « co5 A + ^^sin A. 
A et a ne contenant point A \ donc puisque pour passer des 
fonctions B aux .fonctions E, il suffit de changer co^ A en 

— sinA et mnA en + cos A \ il s'ensuit que si la fonction 

— ou Rco^a.— , par exemple, est représentée par àtf cos A 

A JJ 

+ C sin A , la fonction ^y' . — le sera par jy {—a'sùi A + Ccoi A), 
-et ja sonhne des deux — (B -f f y' E) , le sqn pur 

A 

ou ce qui revi^iit au même y en négligeant les secondes puis- 
sances 4e y' pa^ 

a' eo8(A + }y) + e^sih(A+j/). 

d^QÙ il suit que celte fonction ne serfi autre cfiôse gue la quanr 

^ B C ' ' * • t . " '. 

tité — ou 'Rcosa.r^-y dans laquelle au lieu de A on iqettr<Mt 

A JJ ' 

A+ly'/la partie |j^ on -1(9' — 9) devant être exprimée en 
degrés , minutes et secondes du moyen mouvement du soleil. 

XTS7I. Vn résnltat . semblable aora lieu pour 1^ fonc^tlofi^ 

analogues de B et de £ qui entrent dans nos équations. Soient 
donc 

ce que deviennent les quantités C, F, G, H, lorsqu'on y met 

A + \y fiia pkce de A> ^ nom Mrotu» ^£«9 ^pp^è^ avcni* mî» 

f cos a cos b k la place dç ^ : 

û cas b R 9 9 -^ 

..' I «jJD. jj> ! .1». 






(5a) 

-r- = -T=^ Gi + — ^ Cl 5w a + T-î^ «rt A 
C aD aD .3 

Jjç$ choses, élan t réduites à eet éUt de simplicité ^ oncAlcoIera 
les trois coëfficiens Fs , Ga , H2 , d'après les râleurs 

ï*^ = -^^(«'Fi+yCicosa+yDeosA) 

.99 vl 

Ga = . -J^p (»' Gi +y Ci ««» o +y psinA) (60) 

Ha « 7^7 («' Hi +y Cf tefi^gr b) , •' - - 

■ '.«4 . 

dans lesquelles il feiut se rappeler qu'on a «' =^ ? ( 0' + O » 
y=:l(V — '){et alors on aura simplement , par l'élimination 
de rjDCoimae | , 

n en résulte donc 

m' = M' 4-^ Fa f 

i»' = K + Gap (62} 

P' = Hapî 



• . 



et* en substituant ces valeurs dans l'équation m'*+»'*+p'*= -, 

• r 

elle devient 

- = :^ — i + îip(M'Fa + N'Gfl) + ^*(Pa+G*a+H*2), (53) 

iT* il 

laqablle devra être combinée aved Péquation ordinaire 

r*=R*— p(aRcd^c) + p', 

afin d'en tirer les valeurs de R et p. 

n est visible que ces équations sont de même forme que 
I celles qu'on a obtenues dans le cas de f = ; de sorte que la 

[ solution généjpale , au moins dans l'hypothèse d'une orbite 
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( 33 ) 

parabolique ^ n'offre gaèi*e pluSi.de c<^plioation que. celle qui 
a lieu dans .ce cas . particulier. D'ailleurs les équations précé- 
dentes une fiois résplues^ tout le reste s'achève comme onj'a 
expliqué ci-dessuB, - 

^application à la seconde cqmète de 1781. 

XXXII. Proposons*nons de déterminer Torbite de la seconde 
comète de 1781 , d'après les trois observations suivantes ^ 
réduites pour chaque jour à 8^ sg' 44", temps moyen à Paris. 



9te 



! 



SSB 



SX 



B 



Tnnpê de f^obê. 



Longitude. 



T 



Latitude lu». 



Lieu du 



Log/R. 



NoTémbre. 
A4 



D M S 



D M S 



D M S 



a" 307 

la 3o6 
>' 3o6 



i4 
5i 

4a 



45 
96 



' 



6» 55 
h 39 
b' 3i 



«7 
14 



48 
5s 



a3a 
a37 
i43 



54 
57 



4 



o 4i 



9^94864 1 

9.9944S6 

9.994028 



Ces lieux sont tirés des Mém. de V^cad. des Scienc. année 
1780 1 pag. 67 ; et afin de rendre égaux les intervalles àe temps 
entre les observations, on a interpolé les trois lieux des 19, aa 
et 25 novembre ) ce qui a donné le lieu du a4 tel que nous l'in*- 
serons ici. Ces sortes d'opérations se pratiquent fréquemment 
à l'égard des comètes y et elles sont utiles pour simplifier le 
calcul des élémens. 

Si a"*, a j a sont trois longitudes observées aux temps 
f — m, <)<+/>, respectivement, la longitude 9 pour un tempp 
quelconque t+ z^ intermédiaire entre ceux*là , se calcule par 
la formule 

_ (m+z) (w— jg) ^ 



• = 



'(«~-^)^o^.fOM^^, 



(64) 



Jl en est de mémç des latitudes. 



( 34 ) 

XXXIII. A rinspoction da tabdeaa préoédest ^ on voit que 
le moutremetfit de la comète en longilade est très-pettt ^ cïv^ 
constance peu favorable p<yur rapplication de notre méthode; 
et 8ur-toat pour la détermination du coefficient D; Cependant 
on peut observer que lorsque le moayement en longitude est 
aussi lent .'c|ue dans cet exemple y les observations ne sont pas 
sujettes à de si grandes erreurs , parce qu'alors la comète est 
comparée pendant plusieurs jours à la même étoile, ce qui 
donne des différences assez exactes. 

4 

Pour procéder à rapplîcation de notre méthode / on com- 
la e neera parcaleuler les eoëffieiena G et I> «a mojen ^b»&r* 
mwles (âS) et (ai£r)^ GomBi» il suit>r 

fin 7.83i3893 sùh, 7.4^7670 sin 7.9744880 

to;*g!^.« 9.78016S5 tamgib'*.^ 0.1 Bg^g^Q tangb.. 9.9121887 



.7.6n5646 

Nmb^— o«oo4o884 

— o.ao38do9 

— 0.0079093 



■ ■ ■< 



• 7.5821 Sg9 



• • 



m « 



fien du Qj «37' §7' i" 

180 



7.8866767 

+ 0.0077035 
— 0.0079093 

D =: — o.ooodo6o 



. • A=i= 57* 6/ 4^' ■ ; • 

A — a' = no' 4a' 19' A — o'= m" i4'44' 

I 

sin 9*96g43!27 

iûng B^.... ^- 1^^93999 



51/?. 9.9710027 

ung b*y... 9.78^1653 

* . J '. M" f ' Il ' ■ ■ « I * 

9.7611680 

iVomi... -t- 6.563856 



li > 



• -r- J.345520 

+ 0.563856 
C i= — 0.78^1464 



(55) 

' Où, remarquera que si on eût augnrenté €i d'une miaaRite , It 
eàétàcieni C aaroit subi très-pea de changement; mai^ le 
coë£Bcient D se aeroit réduit à — 0.0000337, valeur neuf fois 
moindre que l'autre. Uite légère augmentation de pins dans la 
longitude a\ auroit fait passer D du négatif au positif; d'où 
l'on voit combien il est à craindre que la valeur de ï) né 
soit pas donnée avec assez de précision par les observations. 
Dana des cas semblables , il conviendra d'éviter tout emploi 
de ce coefficient. Continuons cependant d'appliquer nos for- 
mules, comme si les ébdervationa données étoient exemptes 
é'ertevft. ' . ,• . .^ : 

' -PuisqaeCet D sont â^ ménMiBigney ii «'cncttit qu'où */•< R. 
ti^t^à'dire <tu'iia 19 tioT«tail>i«, époqœ delà seconde obeer. 
vatiou y la tsomète étoit moins éloignée que la ferre du soieiL 
Ensuite pour déterminer r ainsi que p , distance de la comète 
à la terre ^ on résoudra les équations (mg) ;^ mais d'abord il faut 
calculer les coefficiens h et 2Kcos c /d'après lesformnles (fiS)^ 
où Von fera 9 =b= 6** , nyan t soiq d'àjoniter au logarithme de 9 
le log. constant 8. a3558ii« 

A — tf :2= iii*5'$8* C. g>S9Sd^ 

o 

64 

R %.... 9«9g4496 

^ cmù. 9.44S161 .1 :D..*^. 3.686i33 

aR 0.395456 , i : cas b.. o.iiioiS 

(aRco^c)... g. 740617 h i-îSaôôo 

JÙX^tV. Cela posé , en substituant la valettr de R dans letf 
termes tout constans, les éqnatioiîs (219) deviennent 

f = ftf-j i.oSgaôSJ 

1*;;= 0.974658 — f{2Reosc) + f\ 



00s. 9- 556x79 ^fl-^ ^-^69» 

cas ©...-•.. 9.000909 ^ 



(36) 

t)r sachant déjà que r est compris eatmR et Rsm c, on ti'du* 
'Vera aisément par quelques fausses* positiotas : ' ' -^ 

log, r == 9.g9o46a * 

log. p^=± g.7ia736. 

Les vraies valeurs de ces quantités calculées d'après les élé* 
mens connus de Vorbite (vol. cité pag, 71) , sont 

log.r == 9.990071 ^ 

, %-p — 9-709236. . , , . :, 

D'où l'on voit que par ce pœmier calcul la distapce r etst détqr,'{ 
minée y k moins, d'un 1000^ et la distance p à moins d'un 1210*"^^ 
résultat beaucoup plus exact qu'on ne devoit s'y attendre j à 
cause de l'incertitude attachée k la valeur de D , et aussi pareç 
que l'erreur des approximations étant de l'ordre j9% la soluj 
tion ne seaMe devoir être exacte qu'à environ un i4o^| 
lorsque 9 est de 5% qui sont le la^"^ des:58 jours p^is pour unité 

(art;aaan). 

Mais si on exiamine de plus près l'influence de .cette seconde 
cause, on trouvera qu'elle n'est pas à beaucoup près aussi consi-^ 
dérable. En effet, la seconde des équations qu'on vient de 
résoudre , est une équation rigoureuse , qui ne peut être affectée 
que d'une très-légère erreur dans la valeur de cos C. Quant à 
la preqiière équation oii h est la seule quantité qui peut par- 
ticiper a la fois de. l'erreur des observations et de celle de la 
méthode ,• il est aisé de voir , par les valeurs connues de r et 

, , ■ ,' dû 1 dh dr 1 dh 

de 0, qu'on a a tres*peu-pres — ^ = -r.-r-, et — '•= .-7-/ de 

'^ ^ **^ P 10 h r iflo A 

sorte qu'une erreur de j , par exemple , sur A , n'en donne 

qu'une de tî^ sur p, et une de 75^-7 seulement sur r. Cela 

explique suffi^mment, l'exactitude que nous avons obtenue 

dans la solution précédente, sans avoir cependant une valeur 

de D fort approchée. 

XXXY. Venons maintenant aux équations indépendantes 
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* — 

de D et qui doivent donner des réftollats plas certains. H faut 
d'abord calculer les coëfBciens F^ Q, H, d'après les for- 
mules (3i )• Voici le détail de ce calcul : , 



8in {A— a )... 9.9698779 
tang h^ 0.15939:19 



«in (A"^— a*)... 9.^710037 «1» (A — a)... 9*9698779 
$in (A— a* )... 9«9694327 «m(A— a')... 9.9710027 
9 HMg h 0.3182187 tang h' „ 9.780165S 



o.i53654i 

Nomb -f* 1*42447^4 

^. ... 1.7812457 

P = — 0.3567733 



9.7210459 

•— 0.5260729 

— 1.255X728 

1.7812457 



0.0987035 

— X. 2551728 



A— a* = iio* 42' 19* 

A — a'-^iii 14 44 w»(A— a).. 9.9698779 

7 ©o#(A— a').. 9.5691471 
tang b" o*x593929 



• •^ 



2A— û^— a'is 221 57 3^ 



tin,.,....... 9.8250966 

.to^6,..*« .9.9121887 



9.7372853 
— 0.5461 x65 



9.6884179 

+ o»4879979 
-f! 0.1988417 

+ 0.6868396 
— 0.5461165 

* » • 

Qz=: -}- Oii4o723i 



««/i(A— a' ).,. g.9694327 «i>i(A— a*»)... 9.9710027 

tang h 9.9121887 tang 6...... 9.9121887 

tang &»..... 0.1593929 tang h\„^, 9.780x655 



sin{K"^a ).. 9.9698779 
co«(A---^'').. 9.5484643 
tang ^'.«•^., 9.7801653 

• —~— 

9.2985075 

• ^ + 0.1988417 



0.0410143 
+ 1.0990420 



9.6633567 
+ 0.4606548 

+ X. 0990420 



«/i(A— a)... 9.9698779 

tang 5^.. 0,1593929 

ntang &'••. 0.081x953 

0.2 104661 

— 1.6235517 
+ 1.5596768 



. 



4. 1.5596768 H= — 0.0638749 



Il faat ensuite substituer ces valeurs de P/Q^ H, dans T^qua 
lion (3S)*v* • 



*' <• 



(5S) 

Voici le calcul clés coëfliciena dans lesquels la valeur de f sln "¥ 
fi été déduite de Icg. R par la méthode de l'art. XII. 

C 9-892909 P g.SSaSga Q 9.1 48366 

9, 8.93455^ cosi 

z. o ' VA ' " ï.ooiSai . 1.061 531 

I :co^D... 0,111018 90 t.uuij-ii 

8.q3847Q o:6i3gi3 0.209887 

I— .i.* .^ , .«m*.... 8.043643 

cAA A • — n — •'• o-oo55i3 . 

-~- i.o6i5ai "■ ' 8.a5353o 

_ mwnJ. — 4.io3ig4 4- 0.017928 
9.86685i ' —4.16319a 

*-4.i45i64 
^(P+Q* + H') = 20.06817. 

On aura donc l'équation 

^ = o«5i39i7 — /(4.i45a6) + ^•(lo.oSéog) 
r 

log.... 0.617552 1.001478, 

quUl faut combiner avec Féquation déjà trouvée, 

i^Ts 0.974658 — p {^Rcos c) + f* 

log.... 9.7406171 

çt on. trouvera pour solution , 

hg, r ^ 9.990050 
iog. f =! 9.70958a î 

i^aleurs qui ont toute l'exactitude qu'on peut désirer , puisqu'il 
n'y a que très^peu de différence entre ces Valeurs et oellesqu'oi) 
a données ci-dessus d'après les vrais élémens de l'ocbite* 



k 



• — -■-- 
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XXXVI. n faul maintenant déterminer ces élément tels qu'ils 
résultent de lios trois observations : or par les calculs précé-* 
dons on a 



-^P o.6i39i3 

f........ 9.709.58a 


cos h ^ 

-IC Q 

t 


.. 0.909687 
.. 9.70958a 


C09 t 


-H... 9.866851 
f.— 9.70958a 


o.3»^95 


• 


9-9^9469 


• 


/>'..; 9.^76431 


'——^— 2.106177 


fCOth 

♦ C ^ 


— 0.830748 


» 


* 


i>— i 1* 


f lÊfl ?i^.. «.«•.. 


«. o.otiofi7 






..,.-^.— ,.,..... i,^M:àjja 






?*=+ 1.09340A 


— c^4<So5 





Des deux nombres P' et Q' ainsi trouvés , il sera facile, de 
dédoire m' et n' par les formoles (Sa) , qai donnent 

m' =s P' «in A + Q' cos A 

» = — P'cosA + Q'«/ïA. 

ÏInfin les trois autres quantités m^n^p %t ca)cuteront par leè 

formules (^7) > et on aura pour résultat ]çs logarithmes sui- 

vans : 

m g. 881887 m g.GSiSoi 

n .9.715418 ( — /a')....., o.iîi839 

p 9.510753 : p' 9.576433 

où Fon a inAqoé ^Q^ ^' ^^t la seule de ces- six quantités qui 
soit négative. On conclura de-là 

mtil — ib'/i = -— 1.^34955 o.o9i65i 

mp — mp •=^ 0.131674 d'i^d^gg 

,/»/>' — njo = o.6i5j83 9«7&9oo4' 

XXXVn. Fuîsque mv! — m'n est négatif, Te mouvement de 
la comète est rétrograde , et les équations.( 4i ) deviendront : ' 

— ( mp' — m'p) = sin I cos S j/ï n 

— ( np' — np) = sin I sin S j/alT ' 



* 



(4o) 

J}eB deux dernières on tire 

m/?'— Tn'p 9- "94^9 

np' — ri'p 9.789004 

lan^S....... o.66g5o5. 

Ainsi Fangle Sest de 77* 55' 7% ou de aS;^ 55' 7*} or il faut, 
par les équations précédentes , que sin S et cos S soient néga*- 
ti& I puisque sin I ne peut l'être ni |/lin ; donc on a 

S = rkSf 55' 7". 
Connoissant S , on trouvera immédiatement I par la formule 

SI» 9 wn^ 1 =s — —7- -^- , 

qui donne 

I = a6* 59' ^V 5. 

Enfin on aura aussi la distance périhélie par la formule 
y/an 55= • T^^ ^ ^^^ ^^^ **" /flgr. n= 9 , 982474 , ou cette 

distance elle-même 

n = 0.960449. 

L'anomalie déduite de la formule co9^ ? 4 = — 9 sera 

4 = I5•6'46^ 
Avec cette anomalie , on trouve dans la table des oomètes le 

temps correspondant T = lo' • 9697 ; ensuite la formule #=n * T, 
donne #= fo'.3a54. C'est l'intervalle de temps entre l'observa • 
tion du 1 9 novembre et le passage de la comète au périhélie. Pour 
savoir si la seconde époque précède la première ou en est pré- 
cédée , il faut calculer la quantité h , ou seulement en déter- 
miner le signe d'après la formule 

* h = mwl + nri + pp\ 

* 

Or à l'inspection des nombres qui entrent dans cette valeur , 
on reconnoît bientôt que k çst négatif^ et qu'ainsi à l'époque du 



r 



( 4i ) . 

19 BOTembre, & ^9'^!^ temps moyen , la comèie n'aroit pas 
encore atteint son périhélie. Donc en ajoutant à cette date l'in- 
tervalle trouvé lo' 3^54 ou 10'. 7^ 48' 35^, on aura l'instant do 
passage au périhélie : 219 novembre 16^ 18' 19^% ou novembre 

29>6794- ' ■ \^ 

* XXXVm. H reste* a trouver le lieu du périhélie , et d'abord 

il faut trouver la longitude héliocentrique ç de la comète , 

par la formule tangp = — qui donne ^ = 34® 16' 4a'', ou 

p = 9î4» 16' 4a' ; la première est celle qu'il faut choisir , parce 
que sin p et ços f doivent avoir les mêmes signes que n 
et m respectivement , lesquels sont positifs. Enfin la formule 

tang 9 = ^^^ — ^, donne 

^ C03 I 



• ç> 34® 

f — S = i36 



55 7 



ai 



35 



tàng(p^S).... 9.979379 

<^osI 99^9899 

tang<r o.oa948o 

^=i33®3'a5'' 



et par la quantité <r — 4+ S , on trouve ce que les astronomes 
appellent le lieu du périhélie , 

i5® bi' 46*. 
Rassemblant ces divers résultats, on aura les élémens de l'orbite 
de la seconde comète de 1781 , tels qu'ils résultent des trois 
observations données ; nous les mettons ici en comparaison 
avec les élémens corrigés qu'on trouve dans les Mémoires de 
1781, page 71. 



ÉLÉMENS APPROCHÉS. 



,0 



Distance périhélie o . 960449 

Lieu da nœad ascendant 77^ 55' 7 

Inclinaison ^ a6 $9 44 

Lieu du périhélie ]5 5i 46 

Faasige an périhélie : novembre 29' 16^ 18' 19" 



ÉLÉMENS CORRIGÉS. 



o- 960995 
77*» aa* ÔS'' 

16 3 7 

«9* la^" 4** ^6' 



Sens du monrement. Rétrograde, 



f-r. 



•- -■**"^^* ;,.^.,..— ^ 
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XXXIX. La différence de ce^rdeux dyaidinea d^émen» est 
assez petite pour quW ait lieu de s'étonner qu'an premier essai 
de calcul , fondé sur trois observations faites dans riniervalle 
de dix jours seulement ^ ait pu conduire à une connoissanco «i 
approchée de la yéritable orbite. Mais pour comparer encore 
mieux les deux systèmes , nous avons calculé dans chacun 
d'eux les lieux géocentriques de la comète aux époques des 
trois observations. Voici le résultat de ce calcul : 





ÉLÉMBNS APPROCHÉS. 


• 


Époques. 


£A>rt^ 


. calcuiée. 


Long, observée. 


niffér. 


Lotil. caieui. 


Latit, ohêerv. 


Différ, 


i4 

'9 
^4 


D 


M S 


D M S 


S 


D M S 


D M. 


s 


S 


3é7 
3oé 
3o6 


i5 46 
61 37 

4i 59 


3o7 z4 45 
3o6 5i 26 

3o6 42 20 

1 


+ 61 

+ I 
— 21 


55 20 55 

59 '4 49 
3i 3 5o 

• • 1 


55 17 
39 14 
3i 4 


9 

48 

5a 


+226 

+ 1 

— 62 




ÉLÉMENS CORRIGÉ& 

• 


■ 


»9 

24 


3o7 
3o6 
3o6 


17 59 

53 6 

43 28 


3o7 i4 45 
3o(i 5i 26 

* 

3o6 42 20 


+'94 
+100 
+ a8 


* 

55 i5 14 
39 i3 53 
3i 6 24 


55 17 
39 14 
3i 4 


9—115 

48— 55 

62J+ 9» 



/ 



Il résulte de cette comparaison , que nos élémens approehés 
représentent mieux les longitudes que les élémens corrigés , et 
quHls n'ont du désavantage sur ceux-ci que dans la seule lati- 
tude du i4 novembre. Ce résultat inattendu prouve que notre 
méthode a toute l'exactitude nécessaire pour une première 
approximation j mais on devra toujours corriger ^ par des cal- 
culs faits sur des observations éloignées , les élémens donnés 
par des observations peu distantes entre elles j car on ne réussit 
pas toujours à représenter par une même orbite parabolique 
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toatés les observations (Fune comète ; et en cas de dififérerice 
sensible dans les résultats , on doit préférer Forbite qui satis- 
fait le mieax aux observations éloignées. Enfin si la différence 
étoit plus considérable , il faadroit avoir recours à une orbite 
elHptique ou hyperbolique. 

Autre application à la comète de 1769. 

XL. Pour apprécier encore mieux Texactitude de notre 
méthode , nous allons l'appliquer , non à des observations tou- 
, jxrars susceptibles d'erreurs d'où peuvent résulter des compen- 
sations; mais à des lieux de la comète de 1769, calculés d'après 
\qs élémens connus. lies différences qu'on trouvera entre le^ élé- 
mens déduits de ces observations fictives et les vrais élémens , 
seront des erreurs dues tout entières à la méthode , et servi* 
rontàfixerle degré d'approximation auquel elle peut atteindre. 



Epoques. 


1 

Longitude. 


Latitude austr. 


Lieu du 


Log. R. 


Septembre. 
8 à t4^ 


D M S 


D M iS 


D M S 


0.0026648 


a« 101 18 8 


^ as 14 $5 


t66 35 3i 


10 ùi. 


a ria 5i 3i 


b a3 28 i5 


1 68 32 22 


0. 002424 Û' 


la id. 


a' 124 26 47 


V 13 48 3S 


17a 29 20 


0.002I838 



.Voici d'abord le calcul des coëfficiens D et C : 

a^ar = 11^ 33' 23" a — a = 11^ 35' 16' a'~a* ?= aS^S' Sg* 

sin : g.3017509 âin g.SoT^i^d. sin g. 5044435 

tangb'.. 9.6446956 tangb"^. 9.6116907. tangb.. 9-6376970 



8.9464465 
ifqm». -f 0.08839&8 

i + 0.1705480* 



8.9146035 
+ o.o8ai499 



9.â32i4o5 

•— 0.1706636 

+ 0.1705480 

■i ■■ * 

D i= — o.ftaoii55 



^i ifc 
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A :^ 348" 33' aa" 
A — a*= 247*14' 14" A — a = 224^6' 35" 

sin 9*96^7^^9 ^i^ 9.8425oo5 

tangb'.... 9.0446966 tangV*.... 9.6116907 



9.6094805 
JVbmb,.. — 0.4068933 



9.454191a 
+ 0.2845714 



—0.4068933 
+0.21845714 

C=^o.iaa3ai9 



Puisque C et D sont de même signe , il s'ensuit qu'on a r 
comme dans l'exemple déjà apporté j ensuite faisant = 
on aura 

^^ . , ^ /M t . fo.3oio3oo 

COS. 9.7511269 c 9.0875041 

cos b 9.96^4938 R o.oo34a43 

cosc........ 9.7136807 ^ • D 3.9374180 

aR o.3o34549 ^ icosb.. o.o375o6fl 



rf*. 



(aRca5c)... 0.0170749 



h 9.8370467 



XLI. Au moyen de ces valeurs , on aura à résoudre les 
équations 

p = A f — — o . 983394 j 

r^= 1.011226 — pCaRcM c) + p*j 

et comme on sait que r est compris entre R et Visina^ on 
trouvera aisément la solution 

log. r = 9.946165 
log, f = 9.5064^4* 

lies vraies valeurs de ces logarithmes , calculés par les élé- 
mens exacts^ sont 9«946Ô73 et 9.513676^ d'où l'on voit que 
l'erreur est +692 sur le premier logarithme , et — 7263 sur 

le second ^ ce qui répond à une erreur d'un 740^** sur le pre- 

» 



ik^tei^^A 



^ausïà 
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mier nombre , et d'un 6û'"" sur Je second. Si ces premières 
déterminations ne sont pas plus exactes , c'est que 3 A cosc est 
plus grand que R^, mais eu diffère peu 3 de sorte que, suivant 
l'art. XV, les équations à résoudre sont très-près de la limite , 
passé laquelle elles seroient susceptibles de deux solutions. 
Alors les deux couirbéa qui se coupent sont très-près de se 
toucher , et leur intersection Se détermine moins exactement. 
Pour faire cadrer entièrement la solution de nos deux équa- 
tions avec les vraies valeurs de r et de p , il suffiroit de prendre 
ii5o an lieu de ii55 pour la valeur de D, changement qui 
répond A an dixièn>e de seconde sur a^ on sur a\ 

XL II. Il faut passer maintenant au calcul des coëfEciens 
P, Q, H, d'après les formules (3i) : 



fwj (A— a^... 9.9647849 
sin^A — a )... 9.8425oo5 
ikiang b 9.9387270 

9.7460x24 

Nomb + 0.5572017 

— 0.5710873 



P z= -* 0.0138856 



«» (A— a'»).. 9.9647849 «I» (A— a')... 9.84a6oo5 
«m(A — a)... 9.9169326 9in (A — a )... 9.9x69325 
tang b' 9.6446956 tangb^ 9.6116907 



;« » 



9.5264130 
— 0.3360571 

— - o.235o3o2 
— 0.5710873 



9.3711237 
— o.235o3oa 



«»n(A — a ).. 9.9169325 
eo«(A— a').' 9.85625i8 
tang 6^...... 9.6x16907 

9.384875a 
iViojn&*«. — 0.24^59x2, 



Mii(A— a ).. 9.9169325 
co«(A— a®).. 9.58761 74 
tang b' 9*6446956 

9.1492455 

— 0.1410086 
•^ 0,2425912 

— 0.3835998 



À — <»•»=: 249<> 14' 14* 
A — a' 1=224 5 35 

2 A— O**— rt'ss: 1 1 1 19 49 
» « 

ain 9.9691824 

tangb 9.6376970 



9.6068794 

+ 0.4044635 
— 0.3835998 

Qr=: + 0.0208637 
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»« ( A~4s).. . 9.9 1 £93^6 «m(A*^a').,. 9.8426005 «XA— a«)- 9^6^7849 

tan^ b^ S*ii\6^7 U^ b 9-6876970 iang b ••.*.. 9.6376970 

a tang b' 9.9457256 tang b^.... 9.61 16907 «cng 6^ 9-64469^^ 

v^i^MH^^a^i^a^^H^M ^^^^^^^^^^^^^^^^ ^^_^_^_^M_^u.« 

9.4743488 9.091888» 9.247177$ 

î 

Nomb + 0.2980910 T- o.»»356a9 — 0.1766760 

•^ o«3ooâ3JB9 •-" 0.1766760 

H=i — 0.00a z479 — o.3ooa389 

Si l'on veut Toir quel est le réaultat qu^on tireroit des éqoar 
lions (33) et (34)) il &ut multiplier les valeur» d^ PyQ,H« 

R9 
pw -=r- dont le logarithme est 3.4764543^ d'ailleurs par Tar- 

tîcle XII on a log. ( « ^i/i t^ ) = 8. igBtjôiS , ce qui donne î 

P....... ».i4a5ff47 Q g.Sigîgti H 7.332ai4t 

a.4764543 3.4764513 a.4764543 



l-ie* 



R 



0.6I90I90 0.7968457 9.8084684 

.•.9.9975146 isin'^... 8.1969643 



0.61 65336 8.9938100 



Pfbmb,,."^ 4.i35553 + 0.098358 

— 4.135653 



•«-- 4*0371.95 =3 coif. de m. 

On trouve de même | ponr le coefficient de a^^ 

» 

^(WQ^+H*)=i4.i9A4a5^ 

Substituant «s valeurs dojia^ les équations (33) et (34), et 

faisftRt toujours - = o , on aura à résoudre les équations 

a 



r 
1 



-A 



i s 0.494434 — M (2.018598) + •* (7. <x)6a43) 

T 

ô. 851098 



tog.... o.3o6o5o 



r" 



== o . g833g4 + ». 



Or connoissant déjà les limites de r, on trouvera facilement la 

solution 

log. r = g.g4â435 

ïog. » = g.676i48; 

et de la valeur de •« on déduira celle de p = «^--=r • r> 

' "^ D a cos b 

qui donne ibgr* ^ :&=: g • 5i 3ig5. 

Les vraies valeurs de log. r et /^^. ^ étant 9 . 9455^3 et g. 613676, 
on voit que Terreur est — i38 sur le premier logarithme , et 
— 481 sur le second, ce qui ne fait qu'un 3ooo*^ d*erreur 
sur r, et un goo*^ sur p. 

XLIII. Venant ensuite à l'équation (35) ^ on pourra déduire 
ses coëfficiens de ceux de l'équation précédente , a,u moyen de 

a D cos b 
la valeur • = p> p^»/-. j ^^^^ ^^ • ^^j* calculé le coefficient 

pour passer de la valeur de « à celle de p. On aura donc de cette 
manière les équations suivantes en r et f : 

- = o.4g4434 — p(C') + f»(C') 

log.:. o.468oo3 i.i76g34 
r» = i.oiia!^6 — p(aRcMc) + f* 

log.... 0.017076, 
dont la résolution donne 

log. r=3 ^.945619 
log. f =sz g.SiSiaa. 

On voit que l'erreur est encore diminuée sur le premier loga- 
rithme y mais qu'elle est augmentée à-peu-pria d'autant sur le 
second. 
Comme le résultat de ces dernières équations est celui qui 
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mérite le plus de Confiance, lorsqu^on n'a encore aacune con- 
noissance des élémens de l'orbite , noas coptinuerons , d'après 
ce résultat, le calcul nécessaire pour obtenir une première 
approximation vers ces élémens. 

XLIV. Il suffit le plus souvent de connoitre la distance ^ 
périhélie , parce qu'avec cette distance et le rayon vecteur r, i 

on détermine facilement l'anomalie vraie au temps de la seconde { 

observation , et l'instant du passade de la comète par le péri- 
hélie. Or ces premiers élémens étant connus d'une manière 
approchée , on a les moyens de. les rectifier et de déterminer les 
autres élémeps de l'oj^bite , en comparant les observations les 
plus exactes faites à de grands intervalles de temps. 

Si donc on veut se borner à la recherche de la distance péri- 
hélie, il faut faire usage des formules (37) et (4o); savoir, 
n = r— ^A*,et 

k = fcoêb sin (A— a).F+(R— p C05 c) Q' + p sin b.JT j (56) 

formule où l'on a fait pour abréger, 

^ - IÇ-^'^—R ' ^-~¥C~^ 

Q'=ii^Q«.«.^. (57) 

Or dans cet exemple on trouve 

p C05 ô ^w (A— a).... 9.3925:^83 
F..... 7.8983630 

7. 3908913. ... +0.0019538 

R — p cojc....r.... 9.93^7514 
Q' 0.1764585 

0.0982099.... —r. 3537470 
csinb.ïl' 8.2967865^... +0.0198055 

— 1 ,331988=?: it. . 



(49') 

jPui^qiie it est négartif /il A'ensuit qu'au' ïù sejAembrÇ) 'épo<|ue 
de l'observation moyenne y la comète n'avoit pas encore atteint 
son péribélie ; on en déduit 

r 

XLV. Ponr avoir à la fois tous les élémens de l'orbite , tels 
/]u'iJ« peuvent résulter. de ^os trois observations, il faut cal-- 
ciïier les valeurs des six colifiBiciens m^n^p^m'^n^p'^ d'après 
les formules 



m = "R cos A + f cos b cos a 
n s:=z Ti sih A rh p cos b sin a 
p = r^i^ à 



Tri := F sin A + Q' cos A 
n' £=a — P'^^A'+Q'^mA 
p' ^ H' 



et on trouvera leurs logaiithmes comme il suit y en observant 
que ni est la seule de ces quantités qui soit négative : 

m 9^9393289 ( — m').... o.i66a465 

n 8.8788956 ri 9.4lB48i'45 

P g.nSagaS p' 9.i834g37 

De-là on déduit les logarithmes des trois quantités suivantes y 
dont \ts deux premières sont positives et la troisième négative 

(mri^nin) 9-5756968 

{^mp'-^nip) 9.5092008 

— {np'-^rip) 8.4485937 

XLYI. De ce que {mnf-^nin) est positive^ on conclut 
d'abord que le mouvement de la comète est direct ; il feiut donc 
^rendie le sigiie supérieiir dans les' éqiiatighe ( 4i ) , et oti aura 

fnti-^nin î= coèl\/^n. 
mp' — nip =5 sin I cos S j/ain 
np' — n^p = sin I sin S ^/ïn 

rinaj^ction des depx dernières , on voit que sin S doit être 

7 



( 5o ) 

négatif et cos S positif; cle aorte que l'angle S calcalé par la 
formule 

8.4485937 

to/7gS= "^'-"> 9.5093008 

mp — mp .^-_^— 

tangS 8-93939519 

doit être compris entre 270^ et 36o^ ; et comme la tangente , 
qui est négative, appartient également aux deux angles 178^ i' 
45*. 2 et 355' 1' 46''. a, on doit s'arrêter au dernier, 

S = 355^i'45^a. 

Cest le lieu du nœud descendant, suivant la règle de Part. XXV, 
puisque le mouvement est direct et la latitude australe. 
L^inclinaison I se calculera par la formule 

y mp' — m'p 
tang I = r^ — j--f- — 5, 

qui donne I = 4o® 44' li'j. 
Enfin la distance périhélie étant calculée par l'équation 



ti = 



{mn'^^m'ny 

2 C08^ I 



on aura log. n = 9*0913596, ou n s=s o. ia34is6, valeur qui 
s'accorde suffisamment avec celle que nous avons déduite de la 
valeur de t. 
La distance n étant connue, on aura Panomalie 4 par la for* 

mule CM' H = — ^ d'où Ton lire 4 = i36« 4' âo'.4. 

r 

- Cette anomalie répond , dans la table des comètes , au temps 

T = 621'. 749, et ce nombre étant multiplié par n*, on aura 
/= 26'. 9660, qui ajouté au temps de l'observation moyenne 
sept. 10. 5833, donnera ^^p^. 37.5393 ou octobre 7.5393 pour 
rinstant du passage de la comète par le périhélie. 
La longitude héliocentrique f est donnée par la formule 



r 



k 
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tang # = — , où Ton doit observer de rendre le signe de sin f 
m 

conforme à ceIoi.de n , ou celui de cos p à celui de m. On aura 
ainsi f = 4^ 68' 26!' et ^ — S = 9^ 66' 40"^. 8. Ensuite la dis- 
tance ^ de la comète au nœud dont la longitude est S , se calcule 
par la formule 

tang cr = , 

• cos I 

qui donne ^ = i3^ 1' 5a". 8. Enfin la somme S+ <r+4 , puisque 

la comète marche vers son périhélie (art. XXVI) , donnera 

pour le lieu du périhélie, i44^ 8' 8*. 4. 

Voici le résaJcât de tous ces calculs comparé aux vrais élé<- 

mens de Poirbite : 






ÉLÉMENS APPROCHÉS. 



ÉLÉMENS YRAIS. 



Log. dist. périhélie 9.091360 

Passage au périhélie Octobre 7 . 5398 

Lieu du nœud ascendant.. «. . . • 176^ 1' 45^ 

Inclinaison 4^ 44 ^^ 

Lieu du périhélie 144 8 8 






9.090847 

Ociobre.., 7,53io 

1750 y 4o" 

40 47 56 

144 IX 3a 



Sens du mouyement Direct, 



OBSERVATION GÉNÉRALE. 

XliVII. Les deux exemples que nous avons choisis poar 
l'application de notre méthode , sont aussi différens qu'il est 
possible. Dans le premier , le mouvement en longitude est très- 
.lent y et le mouvement en latitude très-prompt. Dans le second, 
au contraire, le mouvement en longitude est très^p rompt, et 
le mouvement en latitude très-lent. Dans l'un , l'intervalle des 
.observations est de 10 jours; dans l'autre, il n'est que de 4 ; 
et cependant le succès de la méthode a été à-peu'-près le même 
dans l'un et dans l'autre. On a obtenu les élémens de l'orbite 
très-peu différens de ceux que donne la combinaison d'un plus 



• --— 
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grand nombre d'obseryations ou d'o^^Y^tionB faites à de plaa 
grands intervalles de temps. 

II y a lieu de prôire que le succès sera, le même dans les 
diverses applications qu'on voudra faire de cette méthode , sauf 
peut-être quelques exceptions qui pourront dépendre ou de 
rinexactitude des observations ou de quelques circonstances 
particulières peu fisivorables , telles qu'une très-petite inclinai- 
son de l'orbite , qui ne permettroit pas d'employer avec sûreté 
les latitudes observées. Le cas le plus défavorable de tous seroit 
celui où la vraie orbite di£féreroit très-sensiblement d'une para^ 
bole y et où en même temps la distance de la comète au soleil y 
lors de la seconde observation , différeroit très-peu de celle dç 
la terre an soleil ; car alors l'emploi des équations (2g) pourra- 
être défectueux à cause de l'incertitude sur la valeur de D, et 
on ne pourroit pas non plus faire usage de l'équation (35), 

puisque - ne seroit pas ass^s petit pour, être négligé.- 
a 

. Au reste, nous avons réduit les formules autant qu'il a été 
possible pour la comm9dité des calculateurs : nous ayons donné 
des règles sûrres pour distinguer , dans le résultat des formules, 
le nœud ascendant du nœud descendant , les époques avant le 
périhélie des époques après le périhélie , et enfin pour prendre 
chaque chose avec le signe qui lui convient, sans avoir recours 
aux figures qui sont cependant utiles pour diriger le calcul. 

Mais on ne sauroit trqp recommander aux calculateurs d'ap- 
porter le plus grand soin à bien déterminer les signes des diffërens 
termes qui entrent dans les formules , et principalement dans 
les valears des .coë£&çiens C , D , P 9 Q ^ H ^ d'après les règles 
connues des sinus, cosinus et tangentes considérés dans les divers 
quarts de la circonférence. Xi'attention , à cet égard ^ est d'autailt 
plus nécessaire , que rien n'y peut suppléer, et que la moindre 
négligence peut faire augmenter considérablement un travail 
déjà long et fastidieux , ou donnera croire mal-à-propos que la 
méthode qu'on a suivie est défectueuse. 



SECONDE PARTIE. 



Méthode pour corriger les élémens de. V orbite connus 

par une première approximation. 

XLVIII. JNous prendrons pour exemple la comète de 1769 
qui a été observée pendant un long espace de temps , avant et 
après son passage par le périhélie, et nous choisirons les trois 
obsefyatious suivantes tirées de la Cométographie de Pingre , 
tome Uy page 368. 



If 



Tempe moyen 
à Paria» 



Longiiude 
de la Comète» 



Août. ii|. 5235a 
SepU. 1 5.61)398 
I Dec; a.2i4i3 



D M S 



39 5S 16 
i4o 39 17 
276 4i 30 



Sa Latitude. 



DM S 



3 17 3o A 
sa 43 34 A 
a3 33 ii5 B 



Lieu du 



Log, R, 



D M S 



\/\2 %\ a6 
173 3i 3o 
260 54 1 a 



• 0052440 
0.0018100 
9 9.J34960 



Supposons que par les premiers essais on a trouvé le l'oga- 
rilhme de la distance périhélie =3= 9*0900000, et l'instant du 
passage au périhélie le 7 octobre à la heures environ; nous 
prendrons pour ces deux élémens corrigés , 

log» n = 9 . 0900000 + • (ï 6000 ) 

Passage au -périhélie octobre 7,5o + t(o*25) 

La correction ^9 ( looop) est exprimée en unités décimales du 
septième orcire ; et la correction t ( ô . 26 ) représente une frac- 
tion de jour qui iroit à ^ de jour ou 6**, si on avoil t = i. Ces 
nombres sont ainsi choisis , parce qu'on suppose que la véri- 
table époque du passage au périhélie, doit être entre 7*35 et 
7.76} et que le vrai log. de la distance périhélie estpareille- 



I 



I 



I 



\ 
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ment compris entre g^oSg et 9,091 ^ ou da moins qu'il s'écdrle 
peu de ces limites. 

XLIX. Gela posé à l'aide de ces deu3t élémeilSj on doit cal- 
culer les anomalies et les rayons recteurs qui ont lieu à Tins-* 
tant de chaque observation. Voici le détail de ce calcul appli- 
qué à la seconde observation : on y verra les précautions à 
prendre pour tenir compte des corrections indéterminées. 

Passage au périhélie... ocl**'* y.So, ou sepl'*~ 3j.5o+r{o.!i5) 
Temps de l'observation sept'*" i5. 69398 



Différence i = 3i'.8o6o2+T(o'.a5) 

Le logarithme de ai .80.602 est i •33S5764 3 pour avoir la partie 
additionnelle due à la correction r (o. âS), on pourroit prendre 
la différence logarithmique de 21806 à 31807^ laquelle est igg, 
et multiplier cette différence par a5o , ce qui donneroit 4g75o ; 
d'où résulte le logarithme entier de < = 1 .3385764 +t (497 5o). 
Mais pour plus d'exactitude, il faut prendre la différence qui 
répond à une unité plus élevée, moitié en dessus , moitié en 
dessous du nombre constant ; par exemple , celle qui a lieu de 
21756 à 2i856. Cette différence est 19916 î on la trouve facile- 
ment, par la disposition des tables, en remontant et descendant 
de cinq degrés dans la même colonne au-dessus et au-dessous 
de 21806. En vertu de cette différence, l'excès logarithmique 
qui répond à t(o. to) sera t (19916) ; donc l'excès qui répond 
à t(o.25) sera r (49790), peu différent de celui qu'on avoit 
trouvé immédiatement On aura ainsi , 

log. / =^ 1 . 3386764 + T ( 49790 ). 
P'ailleurs on a 

I (og* H = 8.635oooo + «( i5ooo). 

Retranchant le second du premier , il restera 

3.7035764 + T (49790 )-^ v(i5ooo). . 
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Cest le log. du temps T employé à parcoarir one égale ano- 
malie 4 dans la parabole dont la distance périhélie est i. 

Le logarithme 2.7035764 répond an nombre 5o5.33a. Pour 
avoir la partie additionnelle , j'observe qae l'unité de diffé-* 
rence entre les nombres 5o4.83 et Co5.83, produit entre.leurs 
logarithmes la différence 8693 ; donc à proportion les diffé- 
rences logarithmiques ^9790 et iBooo , donneront entre les 
nombres les différences 5.795 et i .746. De sorte qu'on aura 

T = 5o5î.33a + t (51 796) — ^ (x' 746). 

Dans la table du mouvement des comètes , on trouve pour 
5o4 jours, l'anomalie i3a* 19' 69% et pour a jours de plus, 
une différence de 4'37'' ou 4^ 4^0. Donc l'anomalie qui répond 
à la valeur de T est 

4=i3a«a3'5r8 + T(ia' 894) — .r (3'885), 

où l'on voit que les corrections sont exprimées en minutes et 
millièmes de minute , ce qui est plus commode pour le calcul 
que des expressions réduites toutes en secondes , attendu qu'il 
faut rapporter les différences à une unité au moins égale à la 
minute, 
n faut chercher maintenant le iogarithme«cosinus de Tangle 
^4 = 66* II' q8'4+ t(6'447) — ^(i'942). 

Ce logarithme, pour la partie connue, est 9.6060431: pour 

avoir l'autre partie , je trouve dans la table qu'une minute de 

plus dans l'angle fait diminuer le Jog.- cosinus de 3865, et 

qu'une minute de moins le fait augmenter de 2863. Le miliea 

de ces différences est 3863 £- pour i' ; donc pour 6^447 la diffé* 

rence sera i846i , et pour i^g^a, elle sera 556 1. On aura 

donc 

log. cos j 4 = 9-6o6o43î — T (18461) + • (556i) 

son double.... 9.3130863 — t (36933) + «(11133) 

log. "n 9.0900000 + «"(loooo) 

Diff. ou log. r = 9.8779138 + T (3693a) — 'ïr (usa) 
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c^est le logaHthtne du rayoti vecteur calcalé d'après la fot«- 

niule (38). 

« 

li. Pour aller plus loin , il est bon de diriger le calcul par 
une figure qui représente à-peu-*près la position relative de la 
terre T, du soleil S , de la comète C, et de sa projection K sur 
. le plan de l'écliplique ( fig. 2 ). 

Longitude du O lyS* 3i' âo* 

Long, de la comète. • . ^ i4o 39 17 

Eloitgalion KTS 3a Sa i3 

Il faut en conclure l'angle CTS entre la comète et le soleil 
(angle qui a été appelé c dans la première partie ) ; on le trouve 
par la formule cos CTS = C05KTS cos CTK , où CTK 

= i = aa^ 43' 34\ 

« 

cos KTS. • • • 9.93422184 
cos b •.**.•. 9*9649015 

cw CTS 9-8891299 CTS =2 39* ly ai'g. 

Dans le triangle CTS, connoissant les deux côtés ST = R, 
es = r^ et l'an des angles opposés CTS , on calculera d'abord 

R 

l'angle TCS par la formule sinTCS = — sin CTS , qui donne 

R.*..«« o«ooi8ioo 
sin CTS.... 9 • B009486 
i:r o.i2ao86a — t(36922) + 4r (112a) 



lA-A 



sinTCS 9.9348448 — T(36gai) 4- «(iiaa) 

De là résultent deux valeurs de TCS (parce que Tringle CTS 
n'est pas ot>tus ) , et il faudra se décider entre l'une ou l'autre , 
d'après la connoissance approchée de l'orbite. Dans ce cas ^ on 
doit prendre la plus grande , savoir 

TCSt= i2aV44'38"8 + t(45' 470) — '^(i' 382); • 
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et comme ôiv a déjà CTS = Sg* i3' ai" g , le troisième angle 
du même trùmgle 

CST=r i8« i' 69' 3— T (45' 470) + «(l'SSa). 
Cela posé , la lalitade héliocentrique C SK = x , se calculera 

«m6«nCST 
par la formule sin a = : — CTS — * 

sinh 9.^869544 

*»»CST.-.. 9.4907547— T (176446) + «(5363) 

i : sin CTS. . . . o. 19905 14' ^^ 

siax... 9.4767605 — T ( 176446) + « (5363) 

De-là résulte x = lo" 54' ^'^ — &c.; mais les corrections de 
l'angle sont innlilés à calculer, et il suffit d'avoir celles du log. 
tosinus. On peut trouver celles-ci, en observant d'après les 
tables , qu'ape différence de 65600 sur le log. sinus , répond à 
une de a434 sur le log. cosinus : on peut aussi multiplier les 
corrections de log. sin x par iang^ x , afin d'en déduire celles de 
log C08 X avec un signé contraire. Ce dernier moyen , qui est le 
plus simple , donnera 

logcos X = 9.9920903 + t(6546) — *» ( ï^S)- 
Maintenant l'angle de commutation KST se calculera par la 

_ ca s C3T 

formule cos KST =5 — : :- • 

*^* cos X 

ca^CST.... 9.9781236+ t(i87Ii) — '^(Seg) 
I : C05 X 0.0079097 — «(6546) +«(198} 

cos KST... 9.9860333 + t(i9i65) — <» (371) 
; Commutât; KST = ,4' :»/ "•3--(37'362)+*(i'i39) 

Longit. de la terre... 355 3i 3o 

Longit béUoc. • = 7«68'4i'3-T(37'36a)+*(i'i39) 

• * * . 

LI Par ces calcula, on a déterminé l'anomalie 4 , la latitude 

8 



/ 
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liéJiocentriqae A ; et la longiluck héliocenlriquè ^ dé la comète 
pour le moment de l'observation du i5 septembre. On ferad» 
semblables calculs pouv les ^pofucs données du. 1 4. août et 
du 2 décembre , et oâ aura les résultati suivons , dans leaqujela 
on a distingué par l'accent *» ceux du i4août^ cl par l'accent ' 
ceux du 2 décembre : 

i4aoâl. j ^^^35715 54.4 + ^ (i«'48o>^«(o'5a3) 

Uoff.sm^''= 8.5267204 •4-'»(i«^4)—«-{g45) 

iSseffembref "^ =?= i33° aa'5e'.8 + *( i»'834)— •(S'SSS) 

pe. luene. ^/^^ ^^ ^ ^ ^ 3767606 - t (i 76446 ) + • ( 5363 } 

» 

2 décembre j 4' =^ 146^ 3/ 1 8* -^ r ( 5'3oo ) — <» ( a'54o) 

i-.!!2L- 1 *'~ '^^ ^ =°-^ "^ "^ C4'7»t> ) + «(0^241 ) 

lep«nta«be. ^/^^_ ^„ ^*_ 9-76ôM79 + t( S^îo ) + « ( »88) 

Lir. CoAceTond- maintenant une aphèce oonceatriqafi au 
soleil , dont la surface soit rencontrée par Je j^n de l'écliptiqtt^ 
et par celui de l'orbite de la comète. Pour mieux saisir la situa- 
tion respective des différens points , on pourra supposer que la 
circonférence de l'écliplique est développée sur un plan (fig. 3), 
et forme la ligne droite Q o O dans laquelle T tepoésâute le 
premier point d'arjV* , rg l'ordfe des signes ^ n le nœud ascen- 
dant, et o le nœud descendant de la comète. A chaque point 
de l'écliplique , comme K, on peut élever une perpendicu- 
laire KC qui représente la latitude de la «omète , en même 
temps que rK est sa longitude j et la suite des points C for- 
mera l'a ligné sinuense'flCoC'n qui représente finferséctiôn 
du plan de l'orbite avec la surface spbérique concentrique âù 
.soleU. Dana cette sorte de projection , les longi^de$ et les lati- 
tudes sont représentées par des grandeurs proporUonnelles ; 
mais les arcs du mouvement héliocentriqœ mnt allérds dans 
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hv^ts proportîous ; ils $ont en général augmenlé» , puisque 
Tare de courbe nCii ne représente qu'âne longueur égale à 08. 
Quoi qu'il en sent , cette figure est très propre a diriger le cal-»- 
cul qui reste à faire pour la détermination de t*orl)ite. 

La distance r^ , supplément -de la longitude du nœud aseen- 
dant û, est d'environ 5"^. Au 14 août, la comète étoît à un 
point C* très - voisin du nœud descendant , puisque la dis- 
fance tK'', complément à 36o* de la longitude ^•, n'est que 
cPenviron S'^.A.u i5 septembre , la comète marchant vers le 
périhélie II, s'est trouvée en C avec la longitude tK d'envi- 
ron 8^ Arrivée au périhélie n ley octobre, elle s'est trouvée 
distante d'environ 3i* du nœud ascendant n. Continuant sa 
route, elle 9. passé le nœud, sa latitude est dévenue boréale, 
et enfin le 3 décembre , elle est parvenue au point C où l'on 
à TK' = 6a^ 3ô' environ , valeur de 36o'— p^SÎ l'orbite delà 
comète est e^cactement parabolique , il ne lui reste à parcourir 
du a décembre 176g jusqu'à l'infini, que l'arc CV, tf^ étant 
égal à nn. Si cette prbite est elliptique , elle continuera son 
mouvementjusqu'à sa prochaine apparition dans la partie o C'^C, 
mais elle ne décrira qu'un angle très-petit dans Un temps con- 
sidérable. 



UII. Après aT<Mr pris celle idée du mouvement héliocen*^ 
trique èé la comète , venons à la détermination tant des ooëflî^ 
ciens t et «- que ftes^ren éiévens de Tarlfitc^* 

Suivant la méthode. ordinaire, on considère le triangle sphé* 
riqXie formé par Tes deux points C*,C^ et par le pôle deJ'éclip- 
tique. Dans ce triangle , on connoîtles'côtés de l'angle au pôle, 
qni sont les e<nnplémetté des làtitàJev &*€% KC ; avec l'ange 
coippris nasuré par V^aec K^'K;*» 36»*— -#*+^ : on peat donc 
déterniaer le troisièioi» côtéC^C. qai doit être égal à la di£^ 
repce des uiMpaUes 4* -^ 4 1 «* ^c^là* réanlte l'équation ,de 
condition 

cw ( J**— +) ^oa$ {*^^f)co$ h CQ$K^+$in h sin a\ (55J 
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^On en obtiendra une semblable de la comparaison des deux 
lieux C et C qui donne 

' cas (4' + 4 ) = cos ( (p — ^' ) cos ^ cos A^— * sin x sin >! j (56) 

• * 

et par le moyen de ces deux équations , on pourra déterminer 
r et ^. 

Cette méthode paroît conduire assez directement.au but j 
mais elle a l'inconvénient de supposer que les trois, points 
O^ C, C sont exactement dans un même plan^ et elle n'offre 
pas les moyens de corriger ou au moins de diminuer Terreur 
qui auroit lieu si cette condition n'étoit pas remplie. 

On n'exprimeroit point la condition essentielle de l'unité 
du plan , en formant par la comparaison des points C et Ç' une 
troisième équation semblable aux équations (55) et (56)«En 
effet , cette troisième équation seroit nécessairement une suite 
des deux autres; car il faudroit que les deux plana C'C% CC 
fissent en G* un angle bien différent de iSo"*, pour que l'arc CC 
mené directement de C en C\ fut sensiblement plus petit que 
la somme des arcs C*C, C**C. Voici donc le moyen qu'on pour- 
roit employer. 

LI Y. On reconnoît que trois points C% C , Cy sont dans le 
plan d'un même grand cercle , si les longitudes de ces pomts 
^% P , p\ et learf» latitudes ^% x, V satisfont à l'équation . 

sin (^— ^**) tang K'+sin (f'^-f) lang9s!'+siH {^^'^')tang x = o , (57) 

ou ce qui revient au mêinei si entre les anomalies 4%49 4^ 
et les latitudes x% x, x', on a l'équation 

• » * * 

• " » 

'«»(4— 4') sin h'+8in (4'^) m ^"+^in (4*-4' )sin?<^ o, (58) 

Si donc ÇU' substitue les valeurs données dans Tarticle LI; 
en observant de prendre 'x' et 4' négativement, on aura une 
troisième équation , laquelle devra être combinée avec les 
équations (55) et (56), afin d'en tirer des valeurs de r 
et de «r , qui y satisfassent aussi bien qu^il est possible. Je dis 
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aussi bien qu^il est possible^ parce qu'ayant trois équaiiond 
l^ur déterminer deux inoonnoes , on ne doit pas prétendre 
satiafaire exaetemenl à ces trois équations , mais seulement 
cherchera rendre leurs premiers membres très-petits. 

LV. Connoissant rei^^ tontes les quantités contenues dans 
le tableau de l'art. LI deviendront entièrement déterminées $ 
on cialculiera ensuite l'ioclinaison I par la formule 

sm (4 — 4 ) 

où Ton aura soin de prendre 4*— 4 (qui représente la dis- 
tance ce*) de même signe que ^ — ^*, parce que rincliuaison 
est toujours supposée moindre que go*. Lorsque I sera déter- 
miné , on aura la longitude du noeud S par la formule 

sin {p — S ) = tang x cos I, (60 ) 

Cette équation, dans laquelle a doit toujours être regardé 
.comme positif y donnera pour l'angle 9 — S deux valeurs com«- 
prises entre o^ et 180^; on en tirera donc deux valeurs de S. 
L'équation sin(^^ — S) =p tang x* cos I , qui a lieu semblable^ 
ment, donnera de même deux valeurs pour S. On aura enfin la 
troisième équation sin ( S— ^') = tangh' cos I ( où l'on a changé 
le signe du premier membre, afin de prendre positivement la 
latitude x% qui est de signe contraire aux deux autres a, a*), 
d'où résulteront encore deux valeurs de S. Or, dans chaque 
système, il est nécessaire que la même valeur de S se retrouve, 
ou exactement ou au moins , à une petite difi*érence près ; celte 
différence venant de ce qu'on n'a pas pu satisfaire exactement 
aux trois équations qui déterminent r et ir. On prendra donc un 
milieu entre ces trois valeurs de S peu. différentes entre elles ^ 
et on aura la longitude d'un nœud auquel on appliquera la 
règle de l'art XXV, pour savoir si c'est le nœud descendant ou 
le nœud ascendant. 

Par les valeurs connue)) de l'anomalie 4 et de la distance 
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périhélie n, on trouvera aisément rinatant du pnMstgb an pértr 
hélie ; le signe positif ou négatif cle ^ —r 9^ indiquera si la miHMr 
yement est direct ou rétrograde ; enfin le lieu du périhélie se 
déterminera par les formules de Part X^VI. Ainsi lacoonoia«' 
sance de tous les élémens de Torbite se déduit aisément de la 
résolution des trois équations (65), (6§) et (S7) ou (66), 
qui donnent les coëfficiens r et «; 

Nous ne faisona qu'indiqua* ceite méthode suis l'appliquer 
à notre exemple, parce que^ quoique sa marche soit assez 
naturelle, cependant elle a rinoonvénient de rie pas réduire 
les erreurs à une même échelle. Deux erreurs égales dans deux 
des. trois équations que l'on considère pourroient répondre à 
des erreurs très-inégales sur lea longitudes ou sur les latitudes ; 
et c'est ce qu'il faut éviter. Nous allons donc exposer une autre 
méthode qui paroit devoir mieux remplir notre but. 

LVI. Supposons que par les formules précédentes on par 
une {Hremière connoissance approchée des élémens de Forbite, 
xm ait trouvé I = 4o'' 5o% S = 355^ 3' (c'est le lieu du nœud 
descendant), et aC ou ^ =n: i6'*48'. 

Appelons pour abréger x la quantité f "— S, on le coténK 
du triangle spliérique rectangle u K^ 9 Q» aura dans ee triangle 
les deux équations 

sin ^ s= êin I ain ^ (61 ) 

t€ag^K:=szcosltaxigw^ {%%) 

iNeus prendrons pour Taienrs corrigées de I et de r les deux 

expressiiMis 

I = 4o* 60' + z (ô'ooo) 

^ = 16*4»' +«(5'ooo); 

el par leur substitution dans l'équation (61) , on aura 

li/il..... 9.8i54854 + js(73oo) 
sin <r 9.4609406 + u (20920) 
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êin ^ 9. 2764310 + z (7300) + « ( aogao). 
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Comparaat oe log. amas à celui de J 'art LI, on a la différeace- 

— 5295+ z (7300) + K(2o9aa) + r (176446)— -^(536^) , 

d'où résulte sur la latitude a Terreur 

E'= — o'5oo + z (l'i i3) + w (3'i89) + r (aB'goo) — ^ (o'8i8). 

Calcoldiit de même le» deux formiitea sinK"" z^z^inlsin ^^^ 
sin k' =i sin I $in 9 , d'après lc& valeurs «•' =» <^ — 4'' + 4> • 
cr' ^= .36o'' — 4' — 4 ^^ ^9 oa trouvera deax valeurs de log. sin ^% 
et log. sin h\ lesquelles comparées aux valeurs semblables déjà 
tJTouvées^ ( art LI ) , donperont sur les angles a*" , ^ lea erreurs 

£"= ô'ïa5 4--j(o'i^)+^»{3îs64) + T(5'657)— *(p'83o) ' 

£'''===— o'ia3 + ;c(4W)—ttC»V&0—''(5'6^)+^(a'i5 

Maintenant le calcul de l'équation $€Êng i^tss.coS'Ï tang <r, donnera 

eosl.: 9.8788748 — z (546o) 

iûngff 9.4798887 + u (22830) 
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tang z^ . . ... 9 . 3587635 ■^-z( S460 ) ■{• u (ai$3o) . 

»=3T52' S'é — z(o'938) + w(3'ga3) 
: a'ailkut» # == 7»58'4ï.3 — T(i7'36a)+«(i'ï39) 

donc-S ou Tu =4* 53' 22^i-z(o'93H)+w(3'9a3)+T(37'36a)-^(i'i 39) 

Les deux équations tang k^'^s:: cas I iang ^"^ et <<90Z^ n!=i ^s ïianga^ 
donneront semblàblement 

Tu5= 4^ 58' 7^5— z<o'i68)+«C3'787H'^(5'33j)---ir(o:46a) 
Tu = 5* 4'3r5 + z(i'96i)+«(6'789>+'r(6'ia8)-^^(7'68o}. 

La seconde vadear de rt> étant retranchée successivement de 
la premièxe'et de la troisième , on aura deux nouvelles différa 
rences ou erreurs , • 

E"= — 4'757~ j8 (o'77o) + r^(o'i36)4-T(4o'e93)-* ^ (o'67o) 
E" = 6'35o+ «(yia9)+wC2'<>02}f t(9'659-)~>(/2ii). 

LVII. Toutes les condiliona d«t problème étant ainsi expri* 
niées, si on vouloit anéantir kAw fois les cinq erreurs £', E"^ 
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E'S E'% E^) on auroit une équation de plus que d'inconnues ^ 
ainsi que le comporte la nature de la question. Mais comme ce 
ne seroit que par un hasard singulier que les valeurs des incon-* 
nues Zj u,r j 9j tirées de quatre de ces équalions, satisferoient 
à la cinquième, et que même dans ce cas, la mélhode que 
nous devons suivre conduiroit au vrai ré^iultat, nous tâcherons 
seulement de diminuer les erreurs de manière à avoir la para^ 
bole qui satisfait le plus exactement qu'il est possible aux trois 
observations données. 

Nous commencerons donc par déterminer r et 'v par la con- 
dition que les erreurs E"' et E'^ soient nulles ; ce qui donnera 

T =s o.i344 + ^(o.o2i43)-+-^(o«ooi3) 
♦ = 1.0606 + jf(o.3278) + i/(6. 2793). 

Substituant ces valeurs dans les trois quantités E% E', Ë'% on 

aura 

E' = :i'.»48 + *(i'.499) + ''(-^'-996) 
E' = o'.oo6.+ «(o'.oSg) +«(3'.o4o) 

E'"= i'.4i8+«(4'.78o)— ttCi'.i67), 

On voit déjà qu'en &isant z et u égaux à zéro y toutes les 
erreurs se réduiront à deuxj l'une E' d'environ 3' 7, l'autre E'" 
de i'4a; erreurs très-tolérables dans la théorie des comètes. 
Mais il est possible de les atténuer encore en cherchant le 
minimum de la somme des quarrés des quantités E', E'", E''\ 

LVIIT; Pour cela, il faut multiplier tous les termes de E' 
par 1*499 coefficient de Zy tous ceux de Ef' par o.oSg, tous 
ceux de Ë^^' par 4.780, et ajouter les trois produits. Opérant 
ensuite de la même manière par rapport aux coefficiens deu, 
on formera les deux équations, 

o =:5 ÎO.148 + «(^5.099) — w(o,^o8) 

0= 5.098 — ^ ( 0.908 ) -i- i/{i9.58o)j 

d'où l'on tire ^ 

» ... JB s= — o.4i44 
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Au tnoy^n Je ces valeurs, les trois erreurs E', E", E"' , 
deviennent 

E' = <i'789, E'= — o'868 , £"'=— o'a56j 

de aorte que la plus grande ne va pas à une minute. 

Substituant ces valeurs de 4; et z^ dans celles de r et de«, 

oif^aura^ > " 

T s= 0.1^40 

«p = 0.8466: 

d'où résulte le log. de la distance périhélie =0.0908466, et 
l'époque du passage au périhélie : octobre 7^53io. 

Des valeurs trouvées pour z et u^ on déduit l'inclinaison 
I = 4o* 47' 55"' 7 , et la distance <r ou oC = i6* 46' 36"!. 
Ajoutant à o C l'anomalie 4 = i3a* sa' 56" 8 + t (ia'894) — 
^ (3' 885) = i3a* ai' i5*4, et retranchant la somme de 1 80% pu 
aura la distance du périhélie au nœud n û = 3o'' Sià' 8' 5. Enfin 
par l'une des expressions de r^ on trouve Tq = 4^^ 56' 19^6; ce 
qui donne la longitude 4u noeud ascendant 176^ 3' 4o^4. 
. Si de cette longitude on retranche la distance n a du péri- 
hélie au nœud , il restera ce que les a^t|:onoipes appellent le 
lieu du périhélie = i44° 1 1' 3i^ g. 

LIX. On peut doiic établir ainsi les élémens de l'orbite de 
la comète dé 1769 , en tant qu'ils résultent des trois observa-* 
tions données : 

Log. dist périhélie 9*0908466 

Lieu du nœnd ascendant • * . 176® 3' ^o^ 

Inclinaison dç l'orbite. 4o 47 56 

Lieu du périhélie i44 1 1 3^ 

Passage au périhélie octçbre 7 . 63xo 

Sens du mouvement direot» 

* ê 

Si d'après ces élémens , on calcule là longitude et la latitude 
géocentriqnes de la comète* k l'instant de chaque observation , 
on trouvera les résultats saivans : 

9 
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Temps moyen 
à Paris. 



Longitude 
calculée^ 



Août. 14.52351 
Sept. 15.69398 
Dec. 2.21 4 13 



39^58' aa' 
140 39 14 
276 4i 22 



Longitude 
observée. 



39*58' 16" 
i4o 39 17 
276 4^ ^o 



Diff' 



>^-^B- 



S7 



+6* 



Latitude 
ûalculée. 



SSB 



3* 16' 44" 

22 44 Si 

23 33 16 



■ I I. I 
Latiif^ 

observée. 



3M7'ïi3' 
S2 43 3I4 
93 33 25 



f * » > 



!>#. 



-89" 

+77 

h- 9 



OÙ l'on voit que les erreurs sont insensibles sur la longitude , 
et ne vont qu'à i'^ sur la latitude. 

Les élémens trouvés par Pingre y d'après les mêmes observa- 
tions ( Cométograpfaie , tome 11, page 38i), représentent bien 
les observations extrêmes du 14 août et du 2 décembre; mais 
sur l'observation du i5 septembre, ils donnent une erreur de 
4-9^ la' sur la longitude, et une de + a' Z']" sur la latitude. 

LX. Au reste, si d'après les élémens que nous ayons trouves, 
ou d'après ceux de Pingre^ on calculoit les lieux Je la comète 
pour d'autres époques où elle a élé observée, on trouveroit 
des erreurs de lâ ou i5 minutes et même plus,.tapt sur la 
longitude que sur la latitude. C'est par cette raison que Pingre 
lui-même trouve différentes orbites en calculant .di^Térentes 
observations j tantôt il trouve l'inclinaison de 4o^ 48' 39'^^ 
tantôt il la trouve de 4o^ 42^ 16", et ainsi des autres ^lémeiis. 

Il n'y a rien à conclure de-là contre l'exactitude des calculs 
d'après lesquels ncrns venons de déterminer l'orbite de la comète 
de 1769; ces calculs -donnent *, entre toutes les paraboles pos- 
sibles , Fune de celles qni satisfont le mieux aux trois obser- 
vations jiomiées. Mais il: est possible*, et il paroit même très- 
probable, d'ftprès ces différences ,' que la vraie orbite de la comète 
de 176g n'est point une parabole; Il faudroit discuter de nou- 
veau lc9 o]i»ser vationa de cette C9mè^ , choiair les plus exacte» 
après les «voir corrigées f de l'aberration et de 4a parallaxe, et 
alors on pourroit essayer de déterminer l'ellipse ou l'hyper* 
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bole qai est fia yraié trajectoire. Le calcul ne seroil pas beau- 
coup plas oompliqoé qaé ka précédens; il y entreroit Félément 

de pla3 -, et troid observations suffiroient toi>jours pour résoudre 

le problème; mai^ il sçroit bon d'en faire entrer au moins 
quatre en comparaison. 

Obêprvaiions sur la méthode précédente. 

LSI. Four avoir une idée nette de la méthode que nous 
proposons 9 et dont nous avons déjtaillé les calculs dans un 
exemple, il faut observer qu'elle est composée de quatre par-' 
tiea distinctes. 

Première partie. Après avoir représenté le logarithme de la 
distance périhélie , et l'instant du passage au périhélie par deuj^ 
expressions contenant chacune une correction indéterminée, 
on circule , d'après ces valeurs , l'anomalie 4 et le rayon veo' 
leur r pour le moment de chaque observation. ' 

Seconde partie. Connoissant tout à-la-fois le rayon vecteui;^ 
.la longitude et la latitude de la comète données par l'observa*- 
tiom , la longitude du soleil «t sa distance à la terre données 
par les tables y 00 calcule la longitude et la latitude héliocen^ 
triques de la comète ; et pour cet effet , il est bon de construire 
une figure qui représente à peu* près les positions respectives 
des différens points , et qui aidera le plus souvent à fixer Findé- 
termination que présente dans certains cas le triangle CTS 
dans lequel on connoît deux côtés , et l'angle opposé à l'un 
id'eax. 

Le xésultat de ces deux premières opérations , dont on peut 
former un petit tableau , donnera , pour l'instant de chaque 
bbservatioti , l'anamalie 4 v 1^ longitude héliocentrique f , et le 
iôg. einue de la Isftitiide héliocentrique a» Il iaut ensuite procé* 
tler à la déteitttttation do plan de l'orbite et de ia position du 
péryiéUe.- 
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Troisième partie. Ayant représenté récliptk|ae par une Ugae 
droite , et l'orbite de la comète y yoe du soleil , par une ligne 
sinueuse; on marquera sur celle-ci Je^ différens lieù^ de la 
comète , d'après les longitudes et latitndes héliocentriqaes déjà 
calculées. On y marquera également le lieu dû périhélie qu'on 
doit déjà connoitre d'une manière approchée par l'une des ano- 
malies. On donnera ensuite des valeurs en partie connues, en 
partie indéterminées, à inclinaison I, et à la distance c- de 
la comète à l'un des nœuds lors de l'observation moyenne. 
Au moyen de ces valeurs , on calculera la latitude a, et la 
quantité k (qui n'est autre chose que la distance 0- réduite à l'éclip- 
tique ) par les formules sin x = sin I sin 0-, tang «== cos I tang ^. 
La comparaison du log. sin x, avec celui qu'on a déjà trouvé 
dans la seconde partie , donne une première erreur £' sur 
l'angle a ; on en trouve une pareille sur les deux autres lati- 
tudes, ce qui fait trois erreurs £', E'^ E". 

Quatrième partie. Des trois valeurs de « correspondantes 
aux trois observations , on déduit trois expressions de la lon- 
gitude du nœud, qui égalées entre elles fournissent deux équa- 
tions entre les quatre coëfficiens inconnus. Tirant de ces équa- 
tions les valeurs de deux des coefficiens , et les substituant dans 
l'expression des trois quantités E', E"^, 'E^*\ il ne restera plus 
qu'à faire en sorte que la somme des quarrés de ces trois quan- 
tités soit un minimum. De-là résultent deux équations qui 
achèvent de déterminer tous les coefficiens , et par suite tous 
les élémens de l'orbite. 

* 

s ^ » 

LXII. On pourra procéder de même , si on veut combiner 
ensemble plus de trois observations ; mais alors il conviendra 
d'établir ddux systèmes d'erreurs. Ijo premier comprendra 
toutes les erreurs K'\, W^Wj &c. relatives ailx .latitudes ; ie 
second comprendra tontes les erreurs ide. la pctoUioA du nœud ^ 
si par les diverses valeurs do « calculées aiix é|[>oques iles diffé- 
rentes observations, on trouve , pour la longitude 4a noeud, 
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les diverses valeurs S', S\ S'\ &c. ; et. qu'on appelle S'+e la 
vraie longitude S du nœud , on aura les erreurs successives 
e' = e\ e" = ^' + S' — S'', e'' = e' + S'— S''', &c. , en nombre 
égal à celui des observations , et e' sera une indéterminée nou- 
velle à joindre aux quatre déjà employées. Dans le système des 
erreurs^', e\ e"', &c., on déterminera e' avec.deux desautiei 
inconnues , de manière que la somme des quarrés de ces erreurs 
soit uri' minimum s on aara ainsi l'expression des trois incon- 
nues en fonction des deux autres, et la substitution de ces 
valeurs étant faite dans les quantités E', Ef\ E''', &c: du pre- 
mier système, il n'y "restera plus que deux inconnues. On 
déterminera enfin ces deux inconnues, parla condition que la 
somme des quarrés des erreurs E' , E" , E*' , &c. soit un mini^ 
mum , et on en déduira les élémens de l'orbite. 

Il seroit plus direct de donner des valeurs en partie connues, 
en partie indéterminées , aux cinq élémens de Torbite , et 
ensuite de calculer pour le moment de chaque observation, la 
longitude et la latitude géocenlriques de la comète. On com- 
pareroit les longitudes et les latitudes calculées avec les lon- 
gitudes etjes latitudes observées, ce qui dônneroit pour chaque 
observation deux différences on erreurs, et on égalerôit ensuite 
à un minimum la somme des quarrés de toutes ces erreurs. Mais 
le calcul conduit de cette manière seroit plus long que celui que 
nous avons indiqué , et il n'en résulteroit qju'un avantage assez 
petit pour la diminution des erreurs, à cause du peu de varia- 
tion que subissent les quantités lorsqu'elles sont voisines du 
minimum. 

Obsewations sur le calcul des corrections indéterminées. 

• 

LXni. Lés tables de logarithmes contiennent tout ce qui 
est nécessaire pour le calcul des corrections indéterminées, 
ainsi qu'on l'a vu dans les exemples précédens. Cependant il 
est bon de remarquer qu'on peut cialculer dans tous les cas les 
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coëfficiens des correclions , sans recoarir aux différences des 
tables qui sont toujours affectées de quelqa'errear. Voici quel- 
ques préceptes à ce sujet. 

Si l'on demande le logarithme du nombre a+^bx, dans 
lequel bx est une correction toujours supposée très-petite par 

m étant le module o« 434^9 , &c. Mais oomme \m logarit^imes 
ont oi^dinairement sept décimales, pour que la correction soit 
exprimée en unités décimales du septième ordre » il faudra, 
aa lieu de m , mettre M = 434^04^ y dont le log. est 6 • 6377Arp^ 

et on aura /og*. {a+ bx) =a + x[ — K 

Réciproquement , ayant log.y s= log. a + x{b)y on en tire 
le nombre y =: a-^ x ( tt). 

Si bx représente un nombre de minutes servant de correc-* 

tion à Tare a; et qu'on* appelle n le nombre constant — 7— , 

{ r' étant le nombre de minutes comprises dans le rayon), dont 
le logarithme est 3. ioi5io4, on aura 

log sin (a-^bx) = log sin a + x{nbcoi-a) 
log cas (a-^bx) = logcoa a^^x (nb tang a) 

logtang(a+bx) = log tanga '\' xi ) 

log cot. (a+bx) = log cot. a — 4C ( — : ) 

formules où les quantités qui multiplient x sont réduites en 
unités décimales du septième ordre , de même que lés loga- 
rilbçies. Réciproquement si on a 

tog wi Y z=z logêin a -^r X {p) 9 il en rétulto y =,a + ' ( — tangaj 

log Ç09 y =? log C05 a -J- jf (r) , il ^n résulte / = a — xi — cot. a i 
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iog- tangy == logiang a -^^ x (c) , il en réavàie jr r=a -f* ' ( — '''^ ^ ^o« ^ ] 
Icg cot. y z= Iog cot. « + * (c) , il en résulle y'=^a — x l — sin a cos. a 1 

la correction ajoutée k la valeur do a étant exprimée en 
minutes. 

Il eat souvent utile de- trouver le Iog. cosinus , lorsqu^on a 
le hg. sinus, ou pice ifersd, sans être obligé de chercher Tare. 
Or par les formules précédentes si l'on a . • 

Iog sin y=log sin a+« (c)*, il en résolté Iog cos y = iog cos a^^xÇc iang* à) 
iogcosyz=:logeosa-^x{fi)^ ilenrétulte logsinyz=zlogéina^^x(ccot.^a)\ 

Avec ces formnles et autres semblables , qu'on peut employer 
suivant les circonstances y on s'habituera aisément au calcul 
des corrections indéterminées qui peut être fortutile dans toutes 
sortes dapproxin^ations. 
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APPENDICE. 

i 

Sur la Méthode dèq moindres quarrés. 

JJams la plupart des questions où il s'agit de tirer des mesures 
données par Tobscrvation , les résnltats les plus exacts qu'elles 
peuvent o£Frir, on est presque toujours conduit à un système- 
d'équations de la forme 

. E = a + ft* + cj^+/c + &c. 

dans lesquelles a,b , c ^ f, &c. sont des coëfiSciens connus , 
qui varient d'une équation à l'autre, et x,y, z, fkc. sont des 
inconnues qu'il faut déterminer par la condition que la Taleiir 
de Ese réduise^ pour chaque équation , à une quantité ou nulle 
ou très-petite. 

Si l'on a aulant d'équations que d'inconnues x,y, z, &c. , 
il n'y a aucune difficulté pour la détermination de ces incon- 
nues , et on peut rendre les erreurs £ absolument nulles. Mais 
le plus souvent, le nombre des équations est supérieur à 
celui des inconnues, et il est impossible d'anéantir toutes les 
erreurs. 

Dans cette circonstance , qui est celle de la plupart des pro- 
blêmes physiques et astronomiques , où l'on cherche à déter* 
miner quelques élémens importans, il entre nécessairement de 
^ "^ l'arbitraire dans la distribution des erreurs • et on ne doit pas 

s'attendre que toutes les hypothèses conduiront exactement aux 
mêmes résultats ; mais il faut sur- tout faire en sorte que les 
erreurs extrêmes , sans avoir égard à leurs signes , soient ren- 
fermées dans les limites les plus étroites qu'il est possible. - 

De tous les principes qu'on peut proposer pour cet objet , 
je pense qu'il n'en est pas de plus général , de plus exact, ni 
d'une application plus facile que celui dont nous avons Sait 
usage dans les recherches précédentes , et qui consiste à rendre 
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minimum la somme ded qaarrés des- erreurs. Par ce moyen , il 
s'établit entre les erreurs une sorte d'équilibre qui empêcliant 
les extrêmes de prévfLloir ^ est très-propre à faire connoître 
l'état du système le plus pro^che de la vérité. 
- La somme desquarrés des erreurs E*-^ E''+ £ *+ &c. étant 

(a + fi« 4- c j' + /■« + &c.)' 
+ (a' + b'x +cy +f'z + &a)* 
+ (a'+ ô"* + c'y +f"z + &c.)' 
+ &c. } 

si l'on cherche son minitnum , en faisant ymer x seule , on 
aura l'équation 

o = /a 6 + xfb* +yfbc-\- zfbf-{- &c. , 

dans laquelle i^^lv fab on entend la somme des produits sem^ 
blables aô + aô'+ a''A^'+ &c. j par /6* la somme des quarréi» 
des coëfficiens de x , savoir &* + 6" + 6''* + &c. , ainsi de suite. 
Le minimum , par rapport ky , donnera semblablemeiit 

o =/ac + xfbc +yfc* + zffc + &c,, 
et le mittinium par rapport à r , 

o =faf+ xfbf + jrfcf+ zff^ + &c., 

où l'on voit que les mêmes toèï^cien^fbcyfbf, &c. sontcom* 
m uns à deux équations, co qui contribue à fEiciliter le calcul. 

£n général y pçur former V équation du minimum par rap^, 
port à l'une, des inconnues , il faut multiplier tous les termes 
de chaque équationproposéfi par le coefficient de l'inconnue dans 
cette équation, pris avec son signe , et faire t^nfi somme de tous 
ces produits» 

On obtiendra de cette manière autant d'équations du mini-- 
mum , qu'il y a d^inconnues , et il faudra résoudre ces équa^ 
tions par les méthodes ordinaire^. Mais on aura soin d'abréger 
tous les calculs, tant des multiplications que de la résolution, 
en n'admettant dans chaque opération que le nombre de chiffres. 

lO 
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entiers ou décimaux que peut exiger le degré d'approximation 
dont la question est susceptible. . 

Si par un hasard singulier, il étoit possible de satisfidre k 
toutes les équationd en rendant toutes les erreurs nulles, on 
obtiendroit également ce résultat par les équations du mini^ 
mum ; car si après avoir trouvé les valeurs de x ^ y^ z ^ &c. 
qui rendent nulles E, E', &a , on fait varier jt, y^ r , &a de 
^XjJ'y^J'Zy 6cc.y il est évident que E* qui étoit zéro deviendra 
par cette variation (afx + bfy + cJ'Zy &c.)*. Il en sera de 
même de E'*, E'* , &c. D'où l'on voit que la somme des quarrés 
des erreurs aura pour variation une quantité du second ordre 
par rapport à J'x^ J'y y &c. ; ce qui s'accorde avec la nature du 
minimum. 

Si après avoir déterminé toutes les inconnues x^jr^ z, &c.^ 
on substitue leurs valeurs dans les équations proposées , on 
connoltra les diverses erreurs E, E% E'^, &a auxquelles ce 
système donne lieu y et qui ne peuvent être réduites sans aug- 
menter la somme de leurs quarrés. Si parmi ces erreurs il s'en 
trouve que l'on juge trop grandes pour être admises, alors on 
rejettera les équations qui ont produit ces erreurs ^ comme 
venant d'expériences trop défectueuses ^ et on déterminera les 
inconnues par le moyen des équations restantes, qui alors don- 
neront des erreurs beaucoup moindres. Et il est à observer qu'on 
ne sera pas obligé alors de recommencer tous les calculs; car 
comme les équations du minimum se forment par l'addition des 
produits faits dans chacune des équations proposées , il suffira 
d'écarter de i'adtlition les produits donnés par les équations qui 
auront conduit à des erreurs trop considérables. 

La règle par laquelle on prend le milieu entre les résultats 
de différentes observations, n'est qu'une conséquence très- 
simfile de notre méthode générale , que nous - appellerons 
Méthode des moindres quarrés. 

En effet , si l'expérience a donné diverses valeurs a'^a^y al^y &c. 
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pour uïie certaine quantité a?, la somme dés quarrés des erreurs 
sera (a'— x)»+ (a"~x)* + (a'''— a:)* + &c. , et en égalant 
cette somme k un minimum , on a ' 
' • o = (a—x) + (a'— 'Jty+ {a'"—x) + &c.; 

j, / . t a + V+a'" + &c. ^^ '- . j^ 

dou résulte jf =* i . ', /i étant le nombre de« 

n 

observations. 

Par^illemeçt^ ai poun déterminer la position d'un point dans 
l'espace 9 on a trouvé, par une première expérience , les coor*^ 
données a% h\ c ; par une seconde , les coordonnées a\ b'\ c\ &c. 
ainsi de suite; soient jp,j^, z^ les véritables coordonnées de ce 
point : alors Terreur de la première e:^périence «era la diatanoet 

du point (a , 6', c' ) au point (x, j^, z)j le quarré de cette 

* ' • ♦ . • • 

distance est • 

et la somme des quarrés semblables étant égalée à un minimum^ 

on en lire trois équations qui donnent spiss— ,j^ = — ,z = — , 

n n n ^ 

n étant le nombre des points donnés par Texpérience. Ces for* 
toules sont les mêmes par lesquelles on trouveront le centre 
de gravité commun de plusieurs masse» égales , situées dans les 
points donnés ; d'où Ton voit que le centre de gravité d'un 
*corps quelconque jouit de cette propriété générab. 

Si on dipise la masse d'un corps en molécules égales et assez 
petkes pour être considérées comme des points j Ta somme dès 
guarr^s des ^istamces des molécules bu Cfntr^ aie grapiié sera 
un minimum. 

On voit donc quie la méthode dça moindires c][uarré$ £ait 
coonoUfe.y en quelque ^ovt» , le oehtne autour duquel viennent 
ae laijger tous les réêullata fouxtua par Texpérience y de manière 
à s'en écarter le moins qu'il est possible. L'ap)plication que nou* 
allons faire de cette méthode à là mesure de la méridienne t 
achèvera de inetire dana tout sen jour sa simplicité et sa 
fécondité. 
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u4pplication à la mesure des degrés du méridien. 

Supposant que le méridien est une ellipse dont les axes sont 
dans le rapport de i à i + <^ ^ si on désigne par D la longueur du 
45*"'* degré, et par S celle de l'arc compris entre les deux lati* 
tudes L et L% on aura par les formules connues ^ et en expri- 
mant h' — L en degrés : 

Q 

S = D (L'—L) — i *D.— «» (L'— L) eo8 (L'+L) ; 
d'où résulte 

S ifin 

+ i«. sin (L'— h) cas (L'+L). 



L— L = 



D 



Comme le 45^^ degré est d'environ sSSoo modules , égaux 

I i+f 

chacun à deux toises , on peut faire -=r = -— - — , C étant une 

D 285oo 

fraction très-petite , et on aura 

S S 270 



L'— L = 



a85oo 



+^.-^^-+«.^^5m(L'— L)co5(L'+L), (a) 



équation qui pour chaque arc dont on connoît la longueur avec 
la latitude de ses extrémités , donnera une relation entre a et C, 
Voici maintenant lés longueurs des différens arcs de la méri- 
dienne de France et les latitudes des parallèles qui les séparent, 
telles qu'elles résultent de l'opération exécutée par les célèbres 
astronomes Delambre et Méchain. 



Lieu . 
de l'observation. 



Sa latitude. 



ArcêcompHtr 
exprimée en modulée^ 



\ 



L' — ti 



■ * ■ m 



L' + L 
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t « 
t 

Panthéon à Paris 



Eyanx. 
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Nous avons donc quatre arcs dont les mesnres étant stibsti* 
tuées successivement dans l'équation (a) ^ fourniront quatre 
équations entre « et C. Mais cooimeices quatre équations ne 
peuvent pas être satisfaites toutes à-k-fois , nous supposerons 
qu'elles ont lieu > en attribuant une certaine erreur à la latitude 
de chaque lieu , et nous appellerons £% E'% &c..les corrections 
additives aux latitudes de Dunkerque , du Panthéon , &c^ Ces 
erreurs n'entrent que dans le premier membre de chaque équa*- 
tion : elles sont trop petites pour affecter le terme multiplié 
par a dans le second membre. Voici donc les équations qui 
résultent des quatre arcs mesurés dans l'opération de la méri- 
dienne , 

E' — E''= 0.0051933 + C (2. 19a) — «(0.563)' 
E''— E'"= o.oo3ioo + C(a.67ii) — «(o.35i) 
E'"—E*^=— 0.001096 +^(2.963)+ «(0.047) 
E'^—E^ = — 0.001808 + C(i.86i) + «(o.fl63) 

Comme il importe de considérer les erreurs séparément , ou 
regardera comme une nouvelle inconnuié l^erreur E'% par 
exemple y et on aura les cinq équations: 

E' = E"+o,oo6oa3 + C(4.864) — «(o.gii) 
E" = E'" + o.oo3ioo+ C(a.67a)— «(o,35i) 
E'''= E^' (*) 

.E'^= E" + 0.001096 — tf (a. 962) — «(0,047) / 

E^ = E*'+ o.oo29Q4-r-C(4.8i3) — «(o.3io) . 

Il faut maintenant faire en sorte que la somme des quarrés de 
ces cinq erreurs soit un minimum, et d'abord cette condition* 
exprimée par rapport à l'inconnue E" dont tous les eoëfficieiia 
sont 1 y donne par l'addition de toutes ces équations : - ' 

o = 5E" + o. oi3i a3 — f (o. aSg ) — « ( 1 .62a) 
donc £"'= —o. 002625 + f (0.648) + «(b^Saé). 

». • ■ 

Substituant celte valeur dans les cinq équations (b)yOn aur^ > 
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E'*fic= 0.003398 + c(4.9t8)-^ «(0.590) 
. B" ss;» 0.000475 + ^(2.726) — * (0.027) 

E*'=s= — o.oo26a5 + C (0.048) + a (0.324 ) (r> 
E*^= <^-*o.ooi529 — ^(d*9i4) + « (0.277) 
E^ === 0.000279^—^(4.765) + ft(o.oi4) 

Pour exprimer ensuite là condition du minimum par rapport 
à f 9 il faut multiplier la première équation par 4^912 coëfiB- 
cient de ^; la seconde^ par 2.720 ; la troisième , par o.o48; 
la quatrième, par —"2-9149 1& cinquième, par — 4.765, et 
égaler à zéro la somme de tous les produits. On opérera sem- 
blablement par rapport à «, et on aura les deux équations : 

0= 0.020983 + f (62.726) — «t(3.83o} ^.^ 

o = ~ 0.003287 — ff( 3.83o) + *(o.53i) ^ ^ 

d'où Ton tire « = 0*00675^ et ^sss 0.0000778 » donc. 

Taplatissement « == -—^r 

i48 

..^ ietle45*-*dfigréD«:iî^ 

• 

L'aplatissement déterminé par la longueur du pendule et par 
quelques phénomènes astronomiques , n'est que de -^^ , et le 
45*^ degré tel qu'on Fa déduit de la comparàisoi;! des mesures 
fiâtes en France avec les mesures faites au Pérou , est de 285o4* lo. 
C'est sur ce dernier résultat qu'est fondée la détermination défi- 
nitive du mètre : il derroit être diminué d'énViroil un 45oo^~, 
si on s'en tenoit aux seules mesures exécutées en France ; maid 
l'aplatissement 7^ tipp peu d'acdord avec celui que l'on oon- 
noit par d'aotrea phénomènes , ne permet pa$ d'adopter oa 
dernier pa^ 

Les valeurs trouvées pour « et f , déterminent l'ellipse qui 
satisfidt aussi exactement qu'il est possible aiix mesures de l'arc 
du méridien compris entre Dunkerque et Barcelonne. Cette 
ellipse est beaucoup plus aplatie que celle qui convient à la 
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figure générale da globe; elle suppose dans les latitudes obser^ 
vées des erreurs que l'on déterminera en substituant les valeurs 
trouvées pour « et f , dans ks expressions de E% Wj &c. : on 
trouvera , en réduisant ces erreurs en secondes , 

La plus grande de toutes ne monte pas à o!\ et la moyenne , 
sans égard. aux signes , n'est que de o''9i. 

Si au lieu de chercher les deux quantités a et C qui con- 
viennent au minimum absolu , on commence par faire la quan- 
tité A égale à l'aplatissement connu j^ , les équations (c) 
deviendront 

E' = o»ooi554 + ^(4'9»î) 
E" = o.oooSgi + C(a.72o) 
E"= — 0.001613 + C(o.o48) 
E"^3s: ~-o.ooo663— C( 2.914) 
E"^ = o.ooo3a3 — ff ( 4.765) 

et on aura pour l'équation du minimum, o^^scoego i o + f (62 . 726)^ 
d'où résulte f = — * o. 0001 436 ; donc 

le 45*- degré = a85oo (1 — = 28504.09 j 

ce qui s'accorde suffisamment avec la détermination adoptée ; 
mais alors les erreurs E', E', &c. exprimées en secondes , 
deviennent 

E' = 3'o6 , E'' = o^oo , E" = — 5' 83, E"= — o'' 88 , E^= 3'62. 

Ces erreurs sont plus fortes qu'elles n'étoient dans le cas du 
minimum absolu ; la plus grande tombe sur la latitude d'Evaux, 
et la moindre , qui est même entièrement nulle , sur celle du 
Panthéon. 

Au reste , les anomalies dans les latitudes , qui sans aucun 
doute ne doivent point être attribuées aux observations , tien- 
nent vraisemblablement à des attractions locales qui agissent 
irrégulièrem^it sur le fil a plomb. Il suffit, pour cela, d'un dé&ut 
d'homoeénéité dans les couches qui avoisinent le point 0ù l'on 
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observe la latitade ; et la même cause qui ra'pptoche le zénith 
apparent du midi on du nord , peut aussi le détourner de quel- 
ques secondes vers Test ou vers Touest ; ce qui explique lés 
inégalités qu'on a aussi observées dans les azimuth^* 

Il résulte de l'existence bien constatée de ces anomalies » que 
la longueur des arcs du méridien est moins propre que celle du 
pendule y à la détermination d'une mesure universelle ; et il 
n'est pas étonnant que des observateurs, d'ailleurs très-exacts , 
ne se soient pas accordés dans les mesures qu'ils ont prises des 
degrés du méridien , puisqu'à raison des attractions locales , les 
latitudes de deux lieux également éloignés de l'équateur ,pour* 
roient différer entre elles de plusieurs secondes. 



Paris , le i& pentâse an i3. 

6 mars j8o5. 
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